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ABSTRACT 


A  theoretical  study  of  metal-semiconductor  (N-type) 
rectifying  contacts  is  developed.  This  study  begins  by 
first  analyzing  previous  models  for  this  type  of  junction. 
Particular  attention  is  given  to  the  Schottky  model  and 
to  the  approximations  it  contains . "  This  model  is  then 
improved  upon  by  taking  into  account  nohuniform  impurity 
ionization  and  the  free  electron  concentration  in  the 
depletion  region.  Using  this  more  exact  model  a 
theoretical  expression  for  the  differential  junction 
capacitance  is  calculated.  The  results  indicate  that 
the  junction  capacitance  as  a  function  of  reverse  bias 
can  be  used  to  accurately  predict  the  doping  concentra¬ 
tion  in  the  semiconductor  material,  but  does  not  yield 
a  correct  measurement  of  the  equilibrium  diffusion 
potential  or  barrier  height. 

The  current  voltage  characteristic  for  this  type 
of  contact  is  also  discussedi.  An  expression  for  the 
I-V  characteristic  of  this  junction  is  derived  based 
upon  a  diffusion  model;  This  expression  is  then 
improved  upon  by  accounting  for  tunneling  and  quantum- 
mechanical  reflection  of  carriers  at  the  junction. 
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CHAPTER  I 

INTRODUCTION 

1.1  Outline  of  this  Study* 

The  purpose  of  this  study  is  to  examine  the  Schottky 
model  for  a  metal- (N-T/pe ^semiconductor  contact,  paying 
particular  attention  to  the  approximations  made  which  limit 
its  accuracy.  The  assumptions  made  which  limit  accuracy  are 
removed  and  a  more  exact  model  proposed,  although  with  great 
sacrifice  in  simplicity*  Poisson* s  equation  based  on  the  more 
exact  expression  for  space  charge  becomes  nonlinear,  but  with 
the  aid  of  numerical  techniques  solutions  are  obtained.  These 
solutions  result  in  more  exact  expressions  for  the  usual 
contact  parameters  than  those  predicted  by  the  Schottky 
model  and  comparisons  between  the  two  theories  are  made. 
Finally,  the  current-voltage  relations  for  the  contact  are 
discussed  and  an  expression  derived  for  the  I-V  characteristic 
of  the  junction  based  on  a  diffusion  model. 

1.2  Historical  Survey. 

The  earliest  systematic  studies  dealing  with  rectifying 
systems  are  generally  attributed  to  Braun.*  Beginning  in 
1874,  Braun  used  a  variety  of  natural  crystals  to  which  he 
applied  base  electrodes  of  various  forms  and  a  point  contact. 


He  then  studied  the  dependence  of  the  total  resistance  of 
the  device  Gn  the  polarity  of  the  applied  voltage  and  on  the 
detailed  surface  conditions  in  the  region  of  the  point  con¬ 
tact.  It  was  also  Braun  who  first  noted  that  the  rectifica¬ 
tion  process  was  located  at  the  contact  itself,  but  he  was 
unable  to  offer  a  general  theory  which  could  predict  his 
findings. 

Possibly  the  first  studies  dealing  with  the  current- 

2 

voltage  character  of  rectifying  systems  were  done  by  Pierce 
in  1910.  He  too  used  natural  crystals  and  studied  the  nature 
of  the  rectification  process  by  an  oscillographic  technique. 
However,  he  was  also  unable  to  offer  a  general  explanation 
for  what  he  had  observed. 

Point  contact  rectifiers  similar  to  those  studied  by 
Braun  and  Pierce  found  wide  application  in  the  early  days 
of  radio  telegraphy  as  detectors,  but  were  not  generally 
understood  and  satisfactory  devices  could  not  be  consistently 
produced.  Although  a  wide  variety  of  rectifying  systems 
were  known  and  the  general  properties  of  each  categorized, 
their  application  always  necessitated  an  adjustment  of  the 
point  contact  ("whisker"  as  it  was  popularly  known)  to  find 
a  sensitive  spot  and  frequent  readjustment  in  order  to 
maintain  useful  rectification. 

The  development  of  a  successful  theory  for  rectifying 

3 

systems  was  the  natural  result  of  Wilson's  interpretation 
of  semiconductlon  based  on  the  band  theory  cf  solids,  which 
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was  first  presented  in  1931*  As  a  first  attempt  at  a  general 

explanation  for  rectification,  it  was  believed  that  tunnel- 

4 

ing  was  responsible  for  the  principal  character.  Wilson 
was  the  first  to  offer  a  quantitative  expression  for  the 
current-voltage  relation  at  the  contact  based  on  tunhel 
theory.  His  results  produced  reasonable  numerical  agree¬ 
ment  with  rectification  ratios  observed  on  cuprous  oxide 
rectifiers;  however,  it  was  later  pointed  out  by  Davyd ov^ 
that  tunnel  theories  predict  a  polarity  of  rectification 
which  is  opposite  to  the  direction  actually  observed.  This 
alone  was  sufficient  cause  to  abandon  the  tunnel  effect  as 
an  explanation  for  rectification,  although  other  discrepan¬ 
cies  also  became  apparent. 

In  1932,  Waibel  and  Schottky^  suggested  that  a  blocking 
layer  of  nearly  stoichiometric  composition  at  the  rectifying 
contact  on  a  cuprous  oxide  rectifier  was  responsible  for  the 
rectifying  character.  In  1938.  Kott^  incorporated  a  form  of 
blocking  layer  along  with  the  effects  of  electronic  diffusion 
and  an  electric  field  to  develop  a  theoreticel  model  which 
could  explain  rectification.  The  Kott  barrier  (See  Appendix 
A),  as  it  became  known,  extended  throughout  the  semi¬ 
conducting  crystal  and  was  to  be  a  special  case  of  a  more 

O 

general  theory  developed  by  Schottky.  In  1939,  Schottky 
suggested  that  the  barrier  associated  with  the  rectification 
phenomena  could  arise  from  stable  space  charges  In  the  semi¬ 
conductor  and  the  presence  of  a  chemically  distinct  layer 
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was  not  necessary  to  explain  its  existence.  Furthermore, 
these  charges  arose  from  the  presence  of  the  metal.  Con¬ 
sistent  with  the  concept  of  stable  space  charges,  Schottky 
devised  a  model  which  could  predict  a  voltage  dependent 
rectifier  capacitance  and  barrier  thickness  and  could  give 
reasonable  agreement  with  experimental  results.  In  1942, 

9 

Schottky  presented  a  final  quantitative  version  of  his 
theory  and  also  assessed  some  of  its  apparent  limitations. 

To  date,  the  Schottky  model  is  generally  accepted  as  the 
approximate  model  for  rectifying  contacts.  Also  in  1942,  a 
model  was  presented  by  Be the, which  was  similar  to  the 
Schottky  version,  but  differed  in  that  it  incorporated  a 
thin  interfacial  layer  of  insulating  material  and  a  slightly 
different  mechanism  for  charge  transport.  (See  Appendix  B 
for  a  complete  analysis  of  the  Bethe  model-) 

The  common  characteristic  of  all  models  developed  thus 
far  was  the  strong  dependence  on  the  difference  of  the 
thermionic  work  function  of  the  metal  and  semiconductor. 
However,  experiments  conducted  on  rectifying  systems  of 
silicon  and  germanium  had  failed  to  show  this  dependence. 

The  apparent  Inconsistency  was  not  explained  until  194 7, 
when  Bardeen"^-  proposed  a  different  mechanism  of  barrier 
formation,  which  was  dependent  on  an  electrical  double  layer 
at  the  free-surface  of  the  semiconductor.  Bardeen’s  theory 
was  able  to  show  that  the  presence  of  such  a  layer  tends 
to  make  the  properties  of  the  contact  independent  of  the  work 


CHAPTER  II 


THEORY  RELATING  TO 
METAL- SEMIC  ONDUC  TOR  SYSTEMS 

2.1  Energy  Band  Diagram  of  a  N-Type  Semiconductor. 

In  order  to  establish  a  common  framework  in  which  a 
discussion  of  the  Schottky  model  will  be  meaningful,  this 
chapter  will  deal  with  the  more  general  theory  relating  to 
metal- semiconductor  systems.  Also,  it  is  of  Importance  to 
note  that  in  this  chapter,  as  well  as  those  which  follow, 
all  discussions  will  be  limited  to  metal«(N-Type )- 
semiconductor  systems  and  any  reference  to  "semiconductor" 
is  meant  to  imply  N-Type  unless  specifically  stated  otherwise 
First,  consider  the  energy  band  diagram  of  an  infinitely 
long  semiconductor  as  sh^wn  in  Fig.  2.1.  This  diagram  shows 
the  relative  positions  of  energy  levels  which  will  be  of 
interest  in  future  discussions.  Here  represents  the 
energy  required  to  transfer  an  electron  at  the  Fermi  level 
of  the  semiconductor  (Efs)  into  free  space.  The  quantity 
is  generally  referred  to  as  the  thermionic  work  function 
and  will  be  a  function  of  the  position  of  the  Fermi  level. 
Since  E^g  is  a  function  of  other  properties  of  the  semi¬ 
conductor  and  will,  in  general,  not  remain  a  constant, 
another  energy  will  be  defined  to  represent  the  energy 


-v'  -  'T'.  1 
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required  to  transfer  s n  electron  from  the  bottom  of  tbs 
conduction  band  (3_)  to  free  space*  ‘Ibis  quantity  is 
called  tee  electron  affinity  and  xill  be  denoted  by  X  . 

i’cit,  consider  the  modification  imposed  cn  fee  band 
structure  dse  to  an  abrupt  termination  of  fee  serf,  conductor 
crystal*  Figure  2*3  shows  the  band  structure  of  a  free 
semi  conductor  surface.  Use  distortion  noted,  or  band  bend- 
ins,  is  due  to  the  presence  of  surface  states  or  allowed 

levels  of  occupancy  which  extend  into  the  energy  gap  {3  }; 

S 

the  levels  being  localized  at  the  surface.  Sine#  the  semi¬ 
conductor  is  N— lype,  the  distortion  in  bend  structure  at 
the  surface  has  created  an  effective  P-Fype  layer,  i*e., 
the  sane  effect  could  have  been  produced  if  acceptor 
inpurities  had  been  In&iffused  from  the  free  surface-  Tims, 
a  P-3J  junction  has  been  formed,  and  due  to  the  polarity  of 
the  surface  potential  (Vs)  the  conduction  band  has  somehow 
been  depleted  of  electrons. 

ifhat  causes  the  surface  potential?  If  ere  examines  the 
physical  situation  at  an  abruptly  terminated  crystal  surface 
the  explanation  becomes  apparent*  Here  the  valence  structure 
of  the  crystal  is  disturbed  since  the  normal  bonding  pattern 
can  no  longer  be  maintained  and  "dangling  bends"  or  unused 

i  p 

valences  are  present  at  the  surface.  Schockley*  has  made 
a  detailed  study  of  how  these  surface  states  arise  free 
atonic  levels  and  under  what  conditions  they  nay  be  expected 
to  appear.  His  general  concision  is  that  one  discrete 


i© 

surface  state,  or  usnsftd  valence ,  can  be  expected  for  each 
surface  a  ten  if  the  lattice  spacing  of  the  original  crystal 
structure  is  sufficiently  snail.  ISsese  stales  fora  a  quasi— 
continuum  over  an  energy  range  within  the  tsmdgap  end  3re 
half  filled  when  the  sarfacs  is  charge-  neutral. 

“Shere  are  additional  factors  which  can  contribute  to 
distortion  in  the  hands  at  the  free  surface.  Cue  could 
imagine  that  sane  of  the  surface  a  tons  could  indeed  be  those 
of  donor  impurities  and  could  thus  donate  electrons  to  empty 
sarface  states.  Still  another  contribution  could  he  due  to 
any  lattice  disturbance  at  the  free  surface,  and  possibly 
the  presence  of  vacancies  or  other  lattice  imperfections. 

“She  important  point  to  be  made  with  regard  to  surface 
states  is  that  they  permit  additional,  levels  of  occupancy 
within  the  handgun  and  give  rise  to  a  potential  harrier  at 
the  free  surface.  Coder  equilibrium  conditions  some  of  the 
electrons  of  the  semiconductor  will  occupy  the  empty  surface 
states  at  the  free  surface,  resulting  in  a  localized  negative 
charge.  Since  sere  electrons  are  non  missing  from  the 
conduction  band  in  the  vicinity  of  the  free  surface,  the 
semiconductor  must  assume  a  net  positive  charge  and  because 
of  the  dimensions  involved,  this  charge  may  be  considered 
as  localized  and  continuous  within  the  region  0<x<5 /.  The 
region  OoxS*  represents  an  inversion  layer  and  is  usually 
referred  to  as  a  depletion  region  since  the  electron  con¬ 
centration  in  this  region  is  essentially  depleted.  In 


addition,  the  positive  cSasge  associated  with  this  region 
causes  the  st rccture  to  become  distorted  as  shown  in 

FIs-*  2,2,  She  potential  denoted  by  ^  Is  ocmsrly  referred 
to  as  the  diffusion  potential  or  equilibrium  harrier  height 
arid  p*(x)  *131  be  need  to  represent  the  negative  potential, 
of  an  electron  at  any  point  within  the  depletion  region 
relative  to  the  ball:  semiconductor. 

Sext.  consider  the  sources  of  the  positive  charge 
associated  with  the  depletion  region: 

(1 }  The  donor  energy  levels®  denoted  fey  3^,  within  the 
depletion  region  are  above  the  Fermi  level  and  are  either 
completely  Ionized  or  at  least  more  Ionized  than  those 
located  In  the  Snalh  semiconductor.  Within  the  depletion 
region  there  will  be  a  net  positive  charge  contributed  by 
each  uncompensated,  Ionized  and  Immobile*  donor  site* 

(2)  Since  the  electron  density  in  the  conduction 
band  decays  rapidly  with  Increasing  energy,  there  will  be 
a  deninlshlng  mincer  of  free-  electrons  as  x  decreases* 
Furthermore,  If  eV„s»KT  then  the  region  C<x<s;  will  be  almost 
uniformly  depleted  of  electrons,  except  In  the  vicinity  of 
x=5J.  Thus ,  one  must  conclude  that  near  x=0  (free  surface) 


♦Impurities  are  decidedly  immobile  in  semiconductor 
crystals  below  a  temperature  of  200-300  °C thus  ionized 
donor  sites  can  be  regarded  as  immobile  below  these  temper¬ 


atures 


the  positive  specs  charge  is  essentially  eli^  while  near 
x=*  the  aval lability  ct  conduction  electrons  Mould  offer  a 
compensating  effect*  For  z?=W  the  net  positive  charge  mist 
of  course  be  zero,  since  in  the  bulk  sen i conductor  no  net 
charge  is  observed. 

{3)  A  final  contribution  of  positive  space  charge  can 
be  attributed  to  the  fact  that  in  the  region  near  the  free 
surface  there  can  be  an  additional  concentration  of  minority 
carriers,  since  the  valence  band  edge  is  relatively  near  the 
Farsi  level.  However,  if  one  considers  a  relatively  wide 
fcsn&gsp  semiconductor  and/or  sufficient  doping  levels,  then 
2j3»HS  and  this  additional  contribution  of  positive  charge 
by  minority  carriers  can  be  neglected* 

The  above  considerations  are  summarized  in  Fig.  2.3. 

The  distributions  shown  are  only  approximate  relationships 
based  on  the  above  discussion,  but  the  concepts  Involved 
will  be  very  useful  in  developing  a  core  exact  mathematical 
model  which  will  be  undertaken  in  a  later  chapter. 

In  Fig.  2.3c  and  Fig.  2.3d  the  depletion  region  has 
been  divided  into  two  approximate  regions  to  emphasize  the 
nature  of  charge  distribution  In  these  areas.  Region  B 
includes  the  area  near  the  surface  and  in  this  region  the 
ionized  donor  atoms  are  the  primary  contributor  to  space 
charge.  Region  A  is  an  area  in  which  the  ionized  donor 
atoms  are  partially  compensated  by  electrons  of  the  conduction 
band.  In  this  region  the  net  charge  density  is  a  function 


(a)  Energy  Diagram 


FIGURE  2,3s  Summary  of  Depletion  Region  Conditions  for  a 

Barrier  Formed  at  the  Free  Surface  of  a  N-Type 


Semiconductor. 
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of  the  distance  x.  Furthermore,  the  barrier  potential  will 
influence  this  distribution  since  Region  B  may  disappear 
entirely  unless  the  barrier  is  sufficiently  high.  By 
sufficiently  high  it  is  meant  that  the  barrier  must  be  high 
enough  to  act  effectively  as  a  barrier  to  conduction 
electrons,  i.e.,  Vd»KT/e. 

2.2  Formation  of  a  Potential  Barrier  at  an  Idealized 

Metal-Semiconductor  Contact. 

When  a  metal  is  brought  into  intimate  contact  with  a 
semiconductor,  a  simple  rectifying  contact  may  be  formed. 

The  condition  necessary  for  rectifying  character  is  that  the 
work  function  of  the  metal  (^m)*  defined  in  the  same  manner 
as  for  a  semiconductor,  exceeds  that  of  the  semiconductor 
(j£s).  This  situation  insures  the  existence  of  a  potential 
barrier  to  electrons  which  in  turn  is  responsible  for  the 
rectifying  nature  of  the  contact.  Also,  the  contact  formed 
may  be  considered  idealized,  since  the  unnecessary  complica¬ 
tions  Introduced  because  of  nonuniform  contact,  surface 
defects  and  the  presence  of  surface  contaminants  will  be 
neglected.  Furthermore,  for  the  present,  the  effects  of 
surface  states  will  also  be  neglected. 

The  formation  of  such  a  simplified,  rectifying  contact 
is  shown  in  Fig.  2.4.  In  Fig.  2.4a  the  energy  band  diagrams 
of  the  two  materials  are  shown  as  they  would  exist  if  both 
materials  were  freshly  cut  and  the  distance  of  separation  (d) 
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(a) 


(b) 


(c) 


(d) 


Before  Contact 
and  Kutuai 
Influence. 


Large  Separation 


Snail  Separation 


Intimate  Contact 


FIGURE  Formation  of  a  Simple,  Rectifying  Metal- 


Seraiconductor  Contact 


Is  large  enough  so  that  any  Interaction  between  the  metal  and 
semiconductor  Is  inhibited.  In  Fig.  2«4b  the  metal  and  seal- 
conductor  are  allowed  to  establish  thermal  equilibrium  with 
a  third  medium  which  is  not  shown,  but  would  be  located 
between  the  two  surfaces*  Cnee  a  continuum  is  formed  and 
thermal  equilibrium  is  established,  the  fact  that  no  net 
current  can  flow  requires  that  the  probability  of  occupancy 
at  any  given  energy  level  must  be  the  same  throughout  the 
continuum;  thus,  the  Fermi  levels  of  the  two  materials  must 
align.*  Furthermore,  if  the  Fermi  level  of  the  semi¬ 
conductor  Is  assumed  to  remain  constant  (an  arbitrary 
reference  point),  then  the  Fermi  level  of  the  metal  must 
rise  relative  to  that  of  the  semiconductor  by  an  amount 
equal  to  the  difference  of  the  two  work  functions,  I.e., 

4-4- 

In  Fig.  2.4c  the  surfaces  are  brought  closer  together. 

As  d  is  decreased  there  will  be  an  Increasing  negative 
charge  built  up  on  the  surface  of  the  metal.  This  buildup 
of  electrons  is  the  result  of  the  conduction  band  of  the 


♦From  the  thermodynamic  point  of  view,  the  Fermi  level 
represents  the  chemical  potential  of  an  electron  in  the  solid, 
Thus,  when  a  continuum  is  formed,  the  chemical  potential  of 
the  mobile  electrons,  and  hence  the  Fermi  level  of  the 
materials  involved,  must  be  the  same. 
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semiconductor  being  at  a  higher  energy  level  than  the  Fermi 
level  of  the  metal,  causing  some  electrons  of  the  semi¬ 
conductor  to  diffuse  into  the  metal*  Furthermore,  the 
electrons  which  leave  the  semiconductor  will  be  forced  to 
reside  on  the  surface  of  the  metal  since  the  metal  is 
already  saturated  with  free  electrons  and  is  an  eaui- 
potential  medium.  Thus,  the  electron  supply  of  the  con¬ 
duction  band  is  depleted  near  the  semiconductor  surface, 
which  in  turn  produces  uncompensated  donor  sites  and  a  net 
positive  space  charge  to  form.  Also,  as  the  difference  in 
charge  between  the  two  materials  is  formed  an  electric  field 
is  produced  which  opposes  further  electron  transfer  to  the 
metal,. 

As  d  is  diminished  further.  Fig.  2,4c.  the  depletion 
region  widens  (W  increases)  and  the  potential  barrier  is 
Increased.  One  should  note  that  will  remain  small  as 
long  as  d  is  fairly  large,  since  most  of  the  potential 
drop  is  across  the  gap  between  the  two  materials.  Finally, 
as  d  approaches  zero  and  intimate  contact  is  made.  Fig*  2.4d. 
an  equilibrium  state  will  be  reached  as  soon  as  the  electric 
field  at  the  contact  prevents  any  further  electron  diffusion 
from  the  semiconductor. 

Cnee  equilibrium  is  established,  a  limiting  value  of 
the  barrier  potential  (V^),  diffusion  potential  (V^),  and 
depletion  region  width  (W)  will  be  reached;  clearly,  the 
limiting  value  of  the  barrier  potential  will  be 


Vb=(#*m"Xe)/e  (2.2.1) 

and  for  the  diffusion  potential, 

Vd=(*fm~^s)/3  (2.2.2) 

Thus,  in  the  absence  of  surface  states,  the  equilibrium 
value  of  the  potential  barrier  for  electrons  in  the  metal 
will  be  the  difference  of  the  metal  work  function  and  the 
electron  affinity  of  the  semiconductor,  whereas  the  height 
of  the  potential  barrier  to  electrons  in  the  semiconductor 
will  be  the  difference  of  the  two  work  functions.  The 
qzzantity  £  -4  is  sometimes  called  the  contact  potential 
difference. 

2.3  Effect  of  Charged  Surface  States  on  the  Contact 
Barrier. 

The  equilibrium  model  for  a  metal- semiconductor  contact 
discussed  in  the  previous  section  has  neglected  the  presence 
of  surface  states  at  the  free  surface  of  the  semiconductor. 
Clearly,  if  these  surface  states  are  present  in  sufficient 
quantity  and  charged  to  a  substantial  degree,  one  must  expect 
some  alteration  in  the  nature  of  the  contact,  at  least  to 
the  extent  of  altering  the  equilibrium  barrier  height  at 


the  contact 


When  the  barrier  is  formed  in  the  same  manner  as  dis¬ 
cussed  earlier  (FI g,  2,4)  an  electric  field  Mill  exist  at 
the  junction  of  the  two  surfaces.  The  establishment  of  this 
electric  field  requires  electric  charge  in  the  two  materials, 
which  must  increase  as  the  distance  of  separation  decreases. 
However,  in  the  presence  of  surface  states  there  is  an 
additional  mechanism  by  which  this  charge  may  be  accommodated 
within  the  semiconductor.  In  addition  to  the  positive  charge 
caused  by  uncompensated  donor  sites,  there  may  also  be  a 
charge  associated  with  the  surface  charge  .induced  in  surface 
states  present  at  the  free  surface  of  the  semiconductor. 

The  density  of  these  surface  states  and  the  amount  of  charge 
they  are  able  to  accommodate  will  determine  the  amount  of 
space  charge  due  to  donor  sites  and  thus  influence  the 
depletion  width.  In  addition,  one  can  no  longer  expect  the 
barrier  height  to  be  given  by  a  simple  difference  of  the 
work  function  and  electron  affinity  of  the  semiconductor; 
Instead,  it  may  well  be  independent  of  these  quantities 
and  depend  solely  on  the  barrier  resulting  from  charged 
surface  states. 

The  influence  of  charged  surface  states  on  the  potential 
barrier  at  a  metal-semiconductor  contact  is  summarized  best 
by  giving  the  conclusions  of  Bardeen* *  in  a  classic  paper 
presented  in  1947,  His  conclusions  will  be  presented  .n  the 
discussion  which  follows. 

If  the  density  of  surface  levels  with  energies  which 
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fall  in  the  fcgndgap  is  sufficiently  high  (greater  than 

4  o  _  2 

1C~^  /ce  e  approximately) ,  there  Mill  fee  an  electrical 
double  layer  formed  at  the  free  surface  of  a  semiconductor. 

This  double  layer  is  formed  from  the  negative  charge  associated 
with  the  charged  surface  states  and  the  positive  space 
charge  associated  with  the  uncompensated  donor  sites  within 
the  semiconductor.  This  double  layer  tends  to  make  the  work 
function  of  the  semiconductor  independent  of  the  height  of 
the  Fermi  level  in  the  bulk  region,  and  thus  independent 
of  the  impurity  concentration. 

The  total  strength  of  the  double  layer  at  a  metal- 
semiconductor  jux-ction  will  be  fixed  by  the  difference  in 
chemical  potentials,  and  thus  dependent  on  the  bulk  prop¬ 
erties  of  the  metal  and  semiconductor.  As  a  consequence, 
the  strength  of  the  double  layer  is  independent  of  the  work 
functions  of  the  material  surfaces  before  contact  is  made. 

The  double  layer  consists  of  the  following  parts: 

(1)  A  double  layer  cf  atomic  dimensions  at  the  metal 
surface • 

(2)  A  double  layer  of  atomic  dimensions  at  the  semi¬ 
conductor  surface . 

(3)  A  double  layer  formed  from  the  surface  charges 
on  the  metal  and  semiconductor,  both  of  atomic  dimensions. 

(4)  A  double  layer  formed  from  the  surface  charge  cf 
atomic  dimensions  and  a  space  charge  extending  to  a  depth 
of  10”6  to  10“4  cm  into  the  semiconductor. 
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The  strengths  of  the  double  layers  nay  be  estimated  as 
follows: 

(a)  If  the  density  of  surface  energy  levels  is  suf¬ 
ficiently  high  (greater  than  10*3  /Cn^t  approximately) ,  the 
double  layer  of  (4)  above  will  be  the  same  as  that  for  the 
free  surface  of  the  semiconductor ,  The  rectification 
properties  will  then  be  largely  independent  of  the  work 
function  of  the  metal,  since  the  difference  in  contact 
potentials  Is  compensated  by  the  double  layer  of  (3)  above* 

(b)  If  the  density  of  surface  energy  levels  is  small 
(less  than  10*3  /Cn2,  approximately),  then  the  double  layer 
of  (3)  above  will  be  small,  and  the  double  layer  of  (*f)  above 
will  be  determined  by  the  difference  in  material  work 
functions# 

(c)  If  the  contact  between  the  metal  and  semiconductor 
is  very  intimate,  it  may  not  be  possible  to  distinguish 
between  the  double  layers  of  (i),  (2),  and  (3)  above#  The 
metal  will  then  tend  to  broaden  the  surface  energy  levels, 
but  If  this  broadening  Is  small  compared  to  the  energy  gap 
of  the  semiconductor,  then  conclusion  (a)  above  will  still 
be  valid. 

(d)  If  the  broadening  of  the  surface  energy  levels  by 
the  metal  is  large,  then  no  conclusions  about  the  space 
charge  of  the  semiconductor  can  be  drawn  from  measurements 

of  contact  potential  differences.  Furthermore,  it  is  possible 
for  all  the  conditions  of  (a),  (b),  (c),  and  (d)  above  to  be 


22 


physically  realizable* 

Since  Bardeen1 s  paper,  there  has  been  extensive  research 

done  in  the  area  of  surface  states,  both  from  an  experimental 

and  theoretical  approach,  in  an  effort  to  be  able  to  predict 

their  influence  on  an  arbitrary  metal-semiconductor  system. 

1 4 

Mead*  has  compiled  a  great  deal  of  the  work  done  on  surface 
states  and  concludes  that  an  arbitrary  metal-semiconductor 
system  may  be  classified  into  two  broad  classes:  (1)  a 
surface  stats  controlled  system,  and  (2)  a  system  in  which 
the  influence  of  surface  states  may  be  neglected.  Further¬ 
more,  his  preliminary  conclusions  indicate  that  the  nature 
of  surface  state  control  for  an  arbitrary  metal  on  an 
arbitrary  semiconductor  may  be  predicted  from  the  nature 
of  the  bonding  mechanism  in  the  semiconductor. 

When  a  metal-semiconductor  system  is  surface  state 
controlled,  a  reasonable  approximation  for  the  barrier  height 
(Eb)  would  be  E  ,  where  Eo  is  the  highest  filled  energy 
level  of  the  surface  states  when  the  free  surface  of  the  semi¬ 
conductor  is  charge  neutral.  Fcr  the  Group  IV  and  III-V 
semiconductors  the  values  of  E  can  be  shown  to  agree  quite 
closely  to  the  relation* ^ 

Ess=Ec-(2/3)Eg  (2.3.D 

or  the  value  of  the  barrier  height  is  approximately  two- 
thirds  of  the  bandgap  energy  (E  ).  (See  Fig*  2,2) 

O 
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when  the  metal-semiconductor  system  behaves  as  If  so 
surface  states  were  present  then  the  value  of  Efe  would  be 
determined  by  the  difference  of  the  metal  work  function  and 
the  electron  affinity  of  the  semiconductor.  In  both  cases 
it  has  been  assumed  that  the  metal- semiconductor  contact 
represents  «n  intimate  contact  between  two  clean  and  uniform 
planar  surfaces, 

Z.k  Influence  of  an  Externally  Applied  Bias  Voltage 
on  a  Contact  Barrier. 

The  discussion  presented  thus  far  has  dealt  entirely 
with  equilibrium  conditions.  If  an  external  voltage  is 
applied  to  the  device  this  equilibrium  is  upset  and  one  would 
expect  the  potential  difference  between  the  two  sides  of  the 
Junction  to  be  Influenced  by  the  external  voltage.  Further¬ 
more,  one  would  expect  the  depletion  region  to  contract  if 
the  potential  has  been  decreased  and  to  expand  If  the 
potential  difference  has  been  increased.  Clearly  then,  the 
equilibrium  model  is  no  longer  applicable  in  the  presence 
of  an  externally  applied  voltage  and  must  be  modified  to 
include  this  influence, 

irfhen  an  external  voltage  is  applied  to  the  device  a 
current  will  flow  and  the  total  voltage  must  be  the  sum  of 
the  contact  potential  and  the  voltages  associated  with  the 
bulk  resistance  and  ohmic  contact  of  the  semiconductor.  In 
most  cases,  at  least  for  reverse  bias  and  low  values  of  \ 
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forward  bias,  the  electric  field  associated  with  the  bulk 
semiconductor  is  small  and  may  be  neglected  along  with  the 
voltage  drop  associated  with  the  ohmic  contact. 

The  following  notation  will  be  adopted  in  discussions 
concerning  bias  and  contact  voltages: 

(1 )  V  will  be  used  to  refer  to  a  bias  voltage  of 

Si 

arbitrary  polarity, 

(2)  V\p  will  be  used  for  Va  when  Va  corresponds  to  a 
forward  biased  condition  of  the  device. 

(3)  will  be  used  for  V  when  V  corresponds  to  a 

i  a  a 

reverse  biased  condition  of  the  device, 

(4)  V.  will  be  used  to  refer  to  that  part  of  the  bias 

voltage  which  appears  across  the  contact  barrier.  Under 
reverse  bias  conditions  and  under  forward  bias  condi¬ 

tions  Vj=Vf-(the  voltage  drops  associated  with  the  bulk 
resistance  and  ohmic  contact  of  the  semiconductor,  when 
these  are  not  negligible). 

Figure  2,5  shows  the  influence  of  Va  on  the  energy 
diagram  of  a  metal-semiconductor  contact.  If  the  Fermi 
level  of  the  metal  is  taken  as  a  reference  and  assumed  to 
remain  constant,  then  eVfi  corresponds  to  the  displacement 
between  the  Fermi  levels  of  the  respective  materials.  It 
should  also  be  noted  that  the  height  of  the  barrier,  as  seen 
by  conduction  electrons  of  the  semiconductor,  increases  when 
Va  is  positive  (relative  to  the  metal)  and  decreases  when 
Va  is  negative.  Since  opposition  to  current  flow  increases 


FIGURE  2.6a;  Energy  Diagram  for  a  Metal-Semiconductor 
Barrier  under  Forward  Bias  Conditions; 
Image  Force  Neglected. 


Bulk  Semiconductor  I  Depletion  Region  |  Metal 
FIGURE  2,5b i  Energy  Diagram  for  a  Metal-Semiconductor 


Barrier  under  Reverse  Bias  Conditions; 
Image  Force  Neglected. 
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when  the  barrier  increases,  a  positive  (with  respect  to 
the  metal)  will  correspond  to  a  reverse  bias  and  a  negative 
Va  to  a  forward  bias.  The  height  of  the  barrier  (E-a)  as  seen 
from  the  metal  is  assumed  to  remain  constant  since  the 
influence  of  any  image  force  (discussed  in  the  next  section) 
has  been  neglected  in  the  diagram. 

The  two  basic  assumptions  on  which  Fig.  2.5  is  based  is 
that  the  charge  within  the  depletion  region  is  not  appreciably 
disturbed  by  the  current  flowing  and  equilibrium  conditions 
are  maintained  within  the  bulk  region  of  the  semiconductor 
and  in  the  metal.  If  this  last  assumption  is  valid  then  the 
Fermi  levels  in  the  bulk  semiconductor  and  in  the  metal  are 
uniquely  defined.  Within  the  depletion  region  the  Fermi 
level  cannot  be  uniquely  defined,  since  when  a  current  is 
flowing  Injection  of  free  carriers  prevents  equilibrium 
from  being  established.  The  Fermi  level  can,  however,  be 
represented  in  the  form  of  a  quasi- Fermi  level  which  repre¬ 
sents  the  electrochemical  potential  for  holes  and  electrons 
separately  as  a  function  of  the  distance  into  the  depletion 
region.  These  quasi-Fermi  levels  can  then  be  used  to  indicate 
a  reference  for  holes  and  electrons  within  the  depletion 
region;  however,  since  the  quasi-levels  are  not  required 
for  the  analysis  herein  they  will  be  omitted  from  Fig.  2,5 
and  no  position  for  the  Fermi  level  is  indicated  within 
the  depletion  region  (0<x<W). 
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2.5  Influence  of  an  Image  Force  on  a  Contact  Barrier. 

In  order  to  demonstrate  the  influence  of  an  image  force, 
consider  an  electron  approaching  a  metal-semiconductor  con¬ 
tact  and  on  the  semiconductor  side  of  the  junction.  This 
electron  will  be  under  the  influence  of  a  potential  which 
exists  inside  the  depletion  region  as  well  as  the  influence 
of  a  Coulomb  attractive  force  as  it  approaches  the  metal 
surface.  This  attractive  force  arises  from  the  presence  of 
an  electron  In  close  proximity  to  the  metal  surface.  It  is 
a  well  known  fact  that  an  electron  of  charge  -e  at  a 
distance  xQ  from  the  metal  surface  will  induce  an  image 
charge  of  +e  at  a  distance  -xQ  inside  the  surface* ^  and 
these  equal  but  opposite  charges  are  then  responsible  for 
the  attractive  force  which  pulls  the  electron  toward  the 
metal  surface.  This  force  is  called,  quite  appropriately, 
an  image  force  and  is  responsible  for  lowering  the  potential 
energy  of  a  conduction  electron  in  the  vicinity  of  the  metal 
surface.  This  use  of  the  conventional  image  force  is  based 

on  the  assumption  that  the  semiconductor  is  acting  as  a 

17 

polarizable  medium  without  free  charge  carriers,  and  seems 
a  valid  assumption  since  the  electron  concentration  near  the 
junction  is  almost  exhausted  of  electrons  for  reasons 
previously  discussed. 

The  reduced  potential  of  an  electron  in  the  vicinity 
of  the  potential  barrier  near  a  metal  surface  requires  an 
alteration  of  the  potential  barrier  near  a  metal-semi- 
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conductor  junction.  Since  the  potential  energy  is  effectively 
reduced ,  the  band  edges  must  bend  down  in  the  vicinity  of 
the  metal  surface  and  eventually  the  potential  of  the 
electron  must  remain  some  finite  and  constant  value  within 
the  metal  -surface.  It  is  important  to  note  that  this 
reduction  in  barrier  height  (eV^)  will  be  independent  of 
the  way  in  which  the  ultimate  barrier  height  (without  image 
effects)  was  determined,  i.e.,  whether  it  was  determined  by 
the  difference  by  surface  states,  or  a  combination  of 

both.  Furthermore,  the  presence  of  image  effects  explains 
the  apparent  discontinuity  of  electron  potential  at  the 
metal  surface,  since  with  the  presence  of  an  image  force  the 
electron  potential  is  no  longer  discontinuous  but  must  decay 
over  some  finite  distance  near  the  junction. 

The  potential  due  to  image  effects  can  be  calculated 
from  the  laws  of  electrostatics,  if  one  assumes  that  the 
semiconductor  acts  as  a  polarizable  medium  free  of  charge 
carriers  within  close  proximity  of  the  metal  surface.  If 
an  electron  of  charge  -e  is  at  a  distance  x  from  the  metal 
surface,  then  by  Coulomb’s  law  the  force  of  attraction  (Fc) 
is 


F  - 
c 


brrS'(2x)a 

s 


l6nS 'xa 

o 


(2.4.1 ) 
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where  e7  denotes  the  high  frequency  value  of  €  •*  The 
s  s 

potential  energy  [e0M(x)]  associated  with  this  force  is  the 
integral  of  the  force  from  the  point  xq  (the  point  at  which 
the  particle  is  located)  to  infinity,  so  that 

f 00 

e&"  (x'j=IF  (x)dx  -  — — -  (2*4,2) 

0  1°  IottS  *x 

o  s  c 

One  immediately  notes  that  Eq.  2.4.2  cannot  be  valid  as 
xo-*0,  since  it  predicts  an  infinite  potential  when  xQ-0 
(see  Fig.  2.6)  and  the  potential  of  an  electron  at.  or  inside 
the  metal  surface  must  have  some  finite  value.  This  apparent 
inconsistency  arises  from  the  approximation  inherent  in 
Eq.  2.4.1#  i.e.,  that  the  distance  of  separation  (xo)  is 
large  compared  to  the  atomic  spacing  of  the  metal  ions.  For 
distances  of  separation  on  the  order  of,  or  less  than  a 
few  atomic  diameters  the  electron  is  most  strongly  influenced 
by  metal  ions  closest  to  it.  The  total  effective  induced 
image  charge  is  still  +e;  however,  the  force  of  attraction 
caused  by  this  induced  charge  is  due  to  many  components, 
each  of  which  is  derived  from  the  neighboring  metal  ions. 


*The  high  frequency  dielectric  constant  must  be  used 
when  dealing  with  image  effects  since  electrons  are  moving 
so  fast  in  the  region  near  the  metal  surface  that  dipole  and 
ionic  polarization  of  the  lattice  does  not  occur. 
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FIGURE  2,6:  The  Potential  of  an  Electron  Due  to  the 

Image  Force  at  a  Metal-Semiconductor  Junction. 


Therefore,  a  more  accurate  model  for  the  force  (Eq«  2.4,1 ) 
is  needed  when  x  becomes  small  and  in  order  to  model  such  a 
force  the  lattice  structure  of  the  metal  must  be  known* 

In  order  to  demonstrate  model  dependence,  consider  a 
simple  cubic  structure  and  an  electron  approaching  the  four 
ions  forming  one  face  along  a  line  passing  through  the  mid¬ 
point  of  the  face  (see  Fig*  2.7)*  Since  the  total  effective 
charge  is  +e»  each  ionic  charge  component  will  be  taken  to 
be  one-fourth  of  the  electron  charge.  The  total  force 

(F  )  will  then  be  the  vector  sum  of  the  force  of  attraction 
c 

resulting  from  each  ion,  or 


e2x 


4rce^(x~+g2/2)3/s 


(2.4.3) 


where  g  is  the  interatomic  spacing  as  shown  in  Fig,  2,7. 
The  potential  energy  of  an  electron  due  to  image  effects 
is  again  found  by  integration,  or 


J 


eg" (x  }=|F„(x)dx  =  - 


4nG^  ( x*!+g s/2 ) 

o  O 


(2.4.4) 


Cne  should  note  that  for  xQ»g/2,  eg"(x0)  will  reduce  to 


(xo»g/2)  (2.4.5) 


which  agrees  with  3q.  2.^.2  except  for  a  factor  of  This 


v/'#*  * 
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factor  of  £  arises  from  the  approximation  that  only  the  four 
nearest-neighbor  ions  exert  an  influence  on  the  approaching 
electron «  Also  of  importance  is  the  fact  that  Ea.  2,4,4  is 
finite  at  xQ=0,  l,e„t 

e3"(0)  =  - - -  (2.4.6) 

(2)3/3Tie^s 

equation  2.4.6  states  that  the  magnitude  of  the  surface 
potential  should  be  inversely  proportional  to  the  inter¬ 
atomic  spacing  in  the  metal.  This  point  is  borne  cut  by 

experimental  findings  for  the  surface  potential  at  a  free 

1 

surface  of  the  alkali  metals*  (Cs,  Hb,  K*  Na,  and  Li)  and 
would  indicate  that  3a.  2.4.6  would  be  the  appropriate, 
although  approximate,  expression  for  the  potential  of  an 
electron  very  near  (x0«g/2)  a  metal- semiconductor  contact. 

Cf  primary  interest  is  the  point  at  which  the  total 
potential  of  the  electron  is  maximum  and  this  maximum  will 
occur  where  the  image  force  on  the  electron  exactly  balances 
the  force  on  it  due  to  the  electric  field  which  exists  in 
the  depletion  region.  Furthermore,  this  maximum  should 
occur  at  a  distance  which  is  greater  than  g,  since  the  major 
decrease  in  electron  potential,  as  predicted  by  Sq.  2.4.4, 
will  occur  within  a  very  few  interatomic  distances  of  the 
metal  surface.  Thus,  one  can  conclude  that  Ea,  2.4.2  should 
be  used  to  approximate  the  upper  limit  for  the  correction 
actually  required  for  image  effects  and  the  corrected  barrier 
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potential  should  appear  as  indicated  in  Pig.  2.8.  The  total 


electron  potential  0(x),  as  shown  in  Fig.  2.8  results  from 
the  addition  of  g*(x)  as  shown  in  Fig.  2.3  and  £"(x)  as 
predicted  by  2q.  2.4.2  for  xi?xmf  or  by  Eq.  2.4.4  as  x-*g. 

The  quantity  is  used  to  denote  the  approximate  correction 
needed  for  image  effects  and  xffi  is  used  to  denote  the  point 
at  which  the  maximum  value  of  0(x)  occurs.  Appropriate 
expressions  for  these  parameters  will  be  deferred  to  a  later 
section  until  an  expression  for  0*(x)  has  been  calculated. 

2.6  Influence  of  Tunneling  and  Quantum-hechanical 

Reflection  on  a  Contact  Barrier. 

The  "tunnel  effect"  arises  from  a  quantum-mechanical 
analysis  of  the  situation  in  which  a  particle  is  Incident 
upon  an  energy  barrier  whose  height  exceeds  the  Kinetic 
energy  of  the  particle.  Classically,  the  particle  would  be 
reflected;  however,  quantua-mechanically  the  particle  has 
a  finite  probability  of  passing  through  the  barrier. 
Quantum-mechanical  reflection  (QKR)  arises  from  a  similar 
analysis  of  the  situation  In  which  a  particle  Is  incident 
upon  an  energy  barrier  whose  height  is  less  than  the  kinetic 
energy  of  the  particle.  Classically,  the  particle  would 
pass  over  the  barrier;  however,  quantum-mechanical  analysis 
predicts  that  there  is  a  finite  probability  that  the  particle 
will  be  reflected. 

Both  of  these  quantum— mechanical  effects  must  be 
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considered  in  the  analysis  of  the  carrier  transport  properties 
of  a  metal-semiconductor  contact.  The  conditions  which 
govern  the  behavior  of  an  incident  electron  on  the  potential 
barrier  associated  with  the  contact  can  be  determined  from 
the  solutions  of  Schroedinger* s  wave  equation.  These 
solutions  would  be  functions  exhibiting  wave  character  and 
would  depend  on  the  electric  potential  in  which  the  electrons 
are  moving.  The  amplitudes  of  the  transmitted  and  reflected 
waves  could  be  inferred  from  the  continuous  nature  of  the 
wave  function  and  its  Instantaneous  spatial  derivative  and 
a  comparison  of  these  amplitudes  would  in  turn  provide  a 
means  of  measuring  tunneling  and  QKR. 

The  usual  procedure  when  considering  the  Influence  of 
quantum-mechanical  effects  on  an  energy  barrier  is  to 
Introduce  the  concept  of  a  quantum-mechanical  transmission 
coefficient  (QKTC),  The  QKTC  for  a  potential  barrier  is 
defined  as  the  ratio  of  the  number  of  electrons  crossing  a 
unit  area  pei  unit  time  in  the  incident  and  transmitted 
waves.  Similarly,  one  could  define  a  quantum-me chan i cal 
reflection  coefficient  (QKRC)  such  that 

QKRC  =  1-  QKTC  (2.6.1) 

and  would  thus  represent  the  ratio  between  the  number  of 
electrons  in  the  incident  and  reflected  waves  on  a  per  unit 
basis. 


Kemble1 y  has  outlined  a  method  for  approximating  the 
QKTC  of  a  parabolic  potential  barrier  in  terms  of  the 
momentum  of  a  particle  incident  from  the  side  corresponding 
to  the  bulk  semiconductor.  The  assumed  form  of  the  barrier 
(parabolic)  could  be  used  to  account  for  the  presence  of 
electron  image  effects  and  would  make  this  method  directly 
applicable  to  estimating  the  QKTC  of  a  metal-* semiconductor 
contact.  The  results  of  Kemble* s  approximation  technique 

are  given  in  Appendix  C. 

20 

Crowell  and  Sze  have  considered  the  problem  of 
calculating  the  QKTC  of  the  potential  barrier  at  a  metal- 
semiconductor  contact  directly  by  using  numerical  techniques 
to  solve  Schroedlnger* s  wave  equation*  The  QKTC  is  calculated 
as  a  function  of  the  carrier  energy  and  effective  mass,  the 
high  frequency  dielectric  constant  of  the  semiconductor  and 
the  shape  of  the  potential  barrier  in  the  vicinity  of  the 
point  where  the  conduction  band  edge  in  the  semiconductor 
merges  into  the  conduction  band  edge  in  the  metal.  The 
mathematical  treatment  involves  numerically  solving  the 
one-dimensional,  time  independent  wave  equation  of  the  form 

p  — 

Ee(x)_ep(x)  li)=  o  (2.6.2) 

where  E  (x)  is  the  electron  kinetic  energy,  m*  is  the 
6  6 

effective  mass  and  eg(x)  is  tho  potential  energy  perturba¬ 
tion  introduced  by  the  barrier.  In  addition  to  assuming  that 
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the  electrons  are  incident  normal  to  the  barrier  (one- 

dimensional  form  of  Eq.  2.6.2)  it  is  assumed  that  m  can 

e 

be  approximated  as  an  average  effective  mass  with  different 
Isotropic  effective  masses  in  the  metal  and  semiconductor. 
Furthermore,  the  image  potential  (g" )  is  assumed  to  have 
the  form  as  given  by  3q.  2.4.2  so  that  the  potential  barrier 
can  be  approximated  by 

p(x)=g*(x)+g,'(x)  (2.6.3) 

for  sufficiently  largo  x. 

A  typical  result  using  the  method  of  Crowell  and  Sze 
is  shown  in  Fig.  2.9.  Their  general  conclusion  is  that  a 
QKTC^0.5  can  be  expected  for  an  electron  incident  on  the 
barrier  with  energy  greater  than  0.05  ev  with  respect  to  the 
barrier  maximum.*  Furthermore,  the  QMTC  Increases  slowly 
with  increasing  energy  and  is  a  rather  strong  function  of 
the  electric  field  at  the  metal** semiconductor  interface. 
Figure  2.9  also  shows  that  tunneling  for  this  particular 
barrier  can  be  appreciable  for  electric  fields  exceeding 
104  v/cm. 

Tunneling  and  QMH  will  be  considered  again  in  Chapter  VI 
when  carrier  transport  across  the  barrier  is  discussed  in 
greater  detail. 


♦Kemble’s  method* ^  predicts  a  QMTC  of  approximately  0.5 
for  electrons  with  the  same  energy  as  the  top  of  the  barrier. 


Electron  Kinetic  Energy  in  Excess  of 
the  Barrier  Maximum  (ev) 

Quantum-Mechanical  Transmission  Coefficient  for 
a  Au-GaAs  (N-Type)  Contact  as  a  Function  of  Electron 
Kinetic  Energy  Relative  to  the  Barrier  Height  for 
Selected  Electric  Field  Strengths  at  tht.  Contact.20 


CHAPTER  III 


THE  SCHOTTKY  MODEL  FOR  A  METAL- (N-TYPE)- 
SEMICONDUCTOR  RECTIFYING  CONTACT 

This  chapter  will  deal  primarily  with  a  mathematical 

o 

treatment  of  the  Schottky  model'  for  a  metal-semiconductor 
contact.  Basically,  the  model  consists  of  assuming  a 
constant  charge  density  in  the  depletion  region  of  the 
semiconductor  and  by  using  Poisson*?  equation  an  estimate  of 
the  electric  field  and  potential  can  be  obtained.  The 
assumptions  on  which  the  model  is  based  will  receive 
particular  attention  as  a  preliminary  step  to  the  develop¬ 
ment  of  a  more  exact  model  which  will  follow  in  the  next 
chapter.  In  addition,  image  effects  will  be  introduced  and 
an  expression  derived  for  the  junction  capacitance  predicted 
by  the  Schottky  model. 

The  energy  band  diagram  for  a  metal-semiconductor 
contact  is  shown  in  Fig.  3«la.  It  is  basically  the  same  as 
Fig.  2.5b  and  is  repeated  here  for  convenience.  It  should 
be  noted  that  the  influence  of  an  image  force  has  been 
neglected  and  also  that  the  junction  should  be  considered 
an  idealized  contact  since  the  effects  of  lattice  imperfec¬ 
tions,  nonuniform  contact  and  the  presence  of  surface 
contaminants  will  be  neglected. 


Contact  under  Reverse  Bias  with  Image  Force. 
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With  reference  to  F3g.  3. la.  the  following  six  state¬ 
ments  describe  the  Schottky  model  for  an  idealised  contact: 

(1)  The  equilibrium  diffusion  potential  (V^)  is  large 
compared  with  KT/e,  i.e,» 

Vd»KT/e  (3.1) 

This  condition  essentially  insures  that  an  effective  barrier 
exists  to  conduction  electrons  of  the  semiconductor.  If  the 
energy  barrier  (eVd)  is  on  the  order  of  ,  or  less  than  the 
thermal  energy  of  an  electron  (KT),  then  the  barrier  is 
easily  surmounted  by  conduction  electrons  and  it  does  not 
act  as  an  ,,effectlve,,  barrier.  Here  "effectiveness"  is 
meant  to  imply  a  measure  of  the  barrier* s  ability  to  prevent 
electrons  of  the  conduction  band  from  crossing  the  junction. 
The  calculations  which  follow  will  be  based  on  the  condition 
that  an  effective  barrier  is  present,  which  necessarily 
limits  application  of  the  model  to  cases  of  small  forward 
bias,  moderate  reverse  bias,  and  the  equilibrium  condition 
stated  in  Eq.  3*1* 

(2)  The  width  of  the  barrier  (W)  is  large  compared  to 
the  wavelength  of  a  conduction  electron.  This  condition 
Insures  an  effective  barrier  against  tunneling,  although 

it  does  not  completely  eliminate  tunneling  from  taking  place. 
Tunneling  must  always  be  considered  as  a  contributing 
factor  to  electron  transfer  at  a  metal- semiconductor  contact, 
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since  by  the  very  shape  of  the  barrier  its  effective  width 
becomes  greatly  reduced  near  the  top. 

However,  if  discussion  is  limited  to  barriers  whose 
width  exceeds  several  thousand  angstroms,  then  tunneling 
may  be  neglected  as  a  major  contributor  to  electron  transfer. 

(3)  The  electron  concentration  on  the  semiconductor 
side  of  the  barrier  is  so  large  compared  to  the  hole  con¬ 
centration  that  it  may  be  considered  a  constant,  i.e., 

n(W)=  n£  (3.2) 


This  allows  to  be  uniquely  specified  and  remain  constant 
for  the  region  x=W. 

(4)  The  concentration  of  uncompensated,  ionized,  and 
immobile  donor  sites  within  the  semiconductor  is  given  by 


(x)  = 


N ,  for  0<x<W 
d 

0  for  W<x<L 


(3.3) 


V. 


and  Nd(W)=nQ,  where  nQ  is  the  equilibrium  concentration  of 
electrons  in  the  semiconductor. 

(5)  The  hole  concentration  at  the  metal-semiconductor 
junction  [p{0)]  does  not  exceed  Nd,  i.e., 


p(0)<Hd+ 


(3.4) 


(6)  The  space  charge  associated  with  the  depletion 


region  is  doriratsd  by  nnccszpensatsd  ionized  donor  sites 


or 


*  f* 

f?(x)dz  «  S^it  »  Mi}2i  (3.5) 

*o 

tender  tbs  above  assertions  egg  easy  T^«n»  Poisson*  s 
equation  to  solve  fer  expressions  for  tbs  electric  field 
snd  potential  distribution  within  tie  depletion  region  of 
tbs  ssElccEdoctoT.  Uns  calculations  wjhjsSa  folios  will 
£ssns  &32S  add  iticnsl  condition  that  fe-yn.  deric®  is  reverse 
biased;  however,  tbs  calculations  sS?na~»d  be  eoually  rslSd. 
for  a  Suia~**l  forward  bias-,  jEne  Uniting  anonnt  of  forward 
bias  will  be  determined  by  the  point  at  wh left  injection  of 
carriers  no  longer  allows  conformance  to  tbs  s<t  statements 
listed  above. 

FtoissGm*s  equation  for  a  planar  contact  *?»»?  xiis  fora 


d2V(z)  -p{x) 

dx-  =  — S^-  l"  j 

s 

in  which  */(x)  is  the  voltage  drop  associated  with  the  contact 
and  p(x)  is  the  space  charge  density  associated  with  the 
depletion  region.  It  is  convenient  to  transform  Eq.  3,6 
into  an  expression  involving  the  negative  potential  of  an 
electron  by  the  relation  which  follows: 


(3.7) 


V(z)  =  -3*(x) 


Tins,  So.  3-o  becomes 


S 


(3-83 


•schere  p*(x)  represents  tie  negative  potential  of  am  electron 
in  tie  depletion  region  witioat  .5 sage  effects.  (See  Fig»  3»i&) 
Since  pfz)  represents  tee  space  charge  density*  it  may 
ce  written  in  general  as 


Pin }  =  e^S^{x5-n{z)-tp(z)] 
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Eowetfer,  for  tie  Schott&y  model  ( statement  #4}  tie  electron 
aid  hole  densities  are  neglected  within  the  depletion  region, 
so  tost  Boissoa’s  equation  becomes 


i25*(z) 

3r2 


(3-iO) 


The  appropriate  boundary  conditions  say  be  easily  seen  from 
Fig«  3>la.  These  are 


(1) 


dg  *  (>: 
dx 


/ 


0  at  x=W 


(3.11) 


(2)  3'(z)  =  0  at  x=W 


(3.12) 


Integration  of  Eq.  3.10  once  with  respect  to  x  yields 


W  in  terms  of  the  barrier  height*  Since 
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£*(  0} 


(3.19) 


It  rollouts  that 


(3.2©) 


JJext,  the  electric  field  (j2(x)3  Is  determines  Is  the 
depletion  region  by  notl ng  that 


a3«(x) 

dx 


-2(3:) 


(3-21) 


Tims,  3(t)  Is  given  by  Ba-  3-15*  or 


2(x)  = 


(O^W)  (3.22) 


Also  of  Interest  is  the  electric  field  at  the  contact  E 

o 

and  fron  Eq.  3-22 

eK,W 

=  ~~  (3.23) 


However,  W  can  also  be  expressed  in  terms  of  the  barrier 
height  (Eq.  3*20)  so  that 


c 

o 


<vv 

O 

s 


r 


(3.24) 
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The  results  of  the  above  calculations  are  susnarlzed  In 


Ksxt,  one  can  determine  the  charge  per  unit  area 
contained  In  the  depletion  region  by  Integrating  the  space 
charge  density  over  the  depletion  width  as  follows: 


J»V 

p(s$dz  =  «3d» 

0 

and  substitution  of  V  fros  So.  3*20  yields 


^  =[&¥s'V 


(3.25) 


(3-26) 


The  depletion  region  evidently  acts  as  a  parallel  plate 
capacitor  since  for  a  seal!  voltage  increase  additional 
charge  will  be  added  near  the  boundary  at  x=W.  Thus,  a 
junction  capacitance  (C  j )  per  unit  area  can  be  defined  as 


(3.27) 


Using  Sq.  3*20,  or  Eq.  3*26,  may  be  expressed  as 
follows:* 


♦For  forward  biased  conditions  (V.+V  )  becomes 

a  r 

(V^-Vj)  and  as  V^— >Vd  the  Junction  capacitance  appears  to 
increase  indefinitely.  This  is 


Bulk  Semiconductor  |  Dsoletira  Region  Metal 

r=¥  x=0 

P IG&BS  3.2a;  Space  Charge  Ds'sity  for  the  Schottky  Model 
of  a  Metal-Semiconductor  Contact . 


FIGURE  >. 2b:  Electric  Field  for  the  Schottky  Model  of  a 
Metal-Seol conductor  Contact. 


FIGURE  3.2c t  Electron  Potential  for  the  Schottky  Model  of 


a  Metal -Semi conductor  Contact;  Image  Force 
Neglected. 
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(3*28) 


In  Chapter  II  the  effect  of  an  loage  force  was  Introduced 
and  the  electron  potential  due  to  the  3 cage  effect  calculated. 
Figure  3.1b  shows  this  influence  and  the  feme  electron 
potential  within  the  depletion  region.  Using  Eq.*s  2.4.2 
and  3*18  the  total  electron  potential,  including  izsage 
effects,  say  be  expressed  as 


g(x)  =  ^(x-*?)3  - 


l6rS*x 
s 


(x>0) 


(3*29) 


Also,  SI  (as  shown  in  Fig.  3.1b)  aay  be  calculated  since 


%  =  e[Vd^Vr-3{xni)]  (3*30) 

where  3(3^)  is  the  maximum  value  of  the  electron  potential 
within  the  depletion  region  when  image  effects  are  included. 
An  expression  for  xm  can  be  found  from  the  fact  that  P(xffi) 
will  occur  at  the  point  where  the  image  force  on  the  electron 


prevented  by  the  fact  that  for  a  small  forward  bias  suf¬ 
ficient  current  flows  to  prevent  from  becoming  comparable 
with  Vd  due  to  the  voltage  drop  associated  with  the  bulk 


semiconductor 


balances  the  force  on  It  due  to  the  electric  field,  or  by 
using  Sq. *s  2.4.1  and  3*22 


eaS-(W-3c  ) 


n 


i6Tiejx^ 
s  n 


(3.3D 


e 


3 


Assuming  that  V»x_  and  solring  for  xE  yields 


a 

“5 


(3.32) 


where  €•  is  the  ratio  of  €*  to  6_.  Using  Eq. *s  3*29,  3.30, 
s  s  s 

and  3*32,  E£  can  be  approximated  as 


*  ”  efe] 


X 

7 


(3.33) 


or  in  terns  of  Dotential  as 


as  6 


re3\'  (V  +7  )' 
w  d'vd+  r^ 


L8na(e,eo  j*e 


S  3 


(3.34) 


by  elininating  EQ  using  Eq0  3*24.  Under  normal  applications 
of  bias  the  lowering  of  the  potential  barrier  to  electrons 
in  the  semiconductor  due  to  image  effects  is  quite  small  and 
is  usually  neglected.  This  will  be  discussed  in  greater  detail 
in  the  next  chapter  and  a  sample  calculation  will  reveal  an 

ft 

estimate  of  the  magnitude  of  lowering  due  to  image  effects. 

The  potential  barrier  as  seen  by  an  electron  on  the 
semiconductor  side  of  the  junction  (Vd  )  may  be  expressed  as 


by  using  Eq.  The  potential  barrier  from  the  metal 

side  of  the  junction  becomes 


and  if  surface  states  are  neglected 


(3.37) 


V.  =  {6  -X  -E»)/e 
bm  me  b 


(3.38) 


The  results  of  this  chapter  can  be  stated  in  the  form 
of  the  following  equations: 


eN 


(1)  e’Cx)  =  2€^(x-W)3 


(volts) 


(2)  p(x)  =  0»(x) - - - 

1 6n€^x 


(3) 


PfsjVV'] 1 


(volts) 


(meters) 
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w 


,  IfWVl5 
1  es  J 


(volts/meter) 


(5)  =  f  2  e  esNd  ( V^4\fr )  ]  ^  (C  oulombs/meter 2 ) 


(*)  C .  = 

j 


re6sKd 


L?(VVr)- 


( farads/meter2 ) 


f  e3N,(V.+y  )"j  * 

(7)  S’  =  e  - A  ■  — 

°  LSira(e*e0)2C  J 


s  s  s 


(electron  volte) 


(8)  Vds  =  Vd  +  Vr  '  Eb/e  (volts) 


(9)  Vta,  =  ‘  WEb)/e 


(volts) 


Comparisons  between  the  Schottky  theory  presented  here 
and  a  more  exact  analysis  will  be  presented  in  Chapter  V. 


CHAPTER  IV 


A  SORE  ACCURATE  MODEL  FOR  THE  IDEALIZED 
METAL- { N-TYFE ) -SEMICONDUCTOR 
RECTIFYING  CONTACT 

The  Schottlcy  model,  as  discussed  in  the  last  chapter, 
is  limited  in  exactness  by  the  assumption  that  all  impurity 
atoms  within  the  depletion  region  are  ionized  and  uniformly 
distributed.  Furthermore,  the  model  assumes  that  all  free 
charge  carriers  are  missing  from  the  depletion  region.  These 
simplifying  assumptions  imply  that  the  total  space  charge 
density  within  the  depletion  region  is  a  constant  and 
outside  this  region  it  is  zero. 

These  restrictions  on  the  space  charge  model  *Till  now 
be  removed  and  a  more  accurate  model  which  includes  both 
incomplete  and  nonuniform  ionization  of  donor  atoms  and  the 
presence  of  free  carriers  within  the  depletion  region  will 
be  substituted.  The  energy  diagram  for  the  metal-semi¬ 
conductor  junction  is  shown  in  Fig.  4.1  and  is  the  same  as 
that  used  for  the  Schottky  mode.  However,  more  quantities 
are  defined  since  they  are  necessary  for  this  more  complete 
derivation  which  will  follow.  One  should  also  note  that  the 
diagram  is  for  reverse  bias  conditions,  just  as  that  for  the 
Schottky  model  and  one  can  again  assume  that  the  expressions 
which  result  will  be  equally  valid  for  small  forward  bias. 


E, — cxl - 
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elV^+Vr) 


T 
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Bulk  Semiconductor  jj  Depletion  Region  q  Metal 

FIGURE  4, Energy  Diagram  for  a  Metal-Semiconductor 
Contact  under  Reverse  Bias;  Image  Force 
Neglected. 
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In  addition,  it  is  assumed  that  any  free  carriers  present  in 
the  depletion  region  are  due  to  the  electrons  furnished  by 
ionized  donor  atoms  rather  than  injected  from  the  metal. 

Then  from  the  standpoint  of  an  idealized  junction  the  model 
should  be  exact,  although  exactness  must  be  relative  since 
tunneling,  surface  effects,  and  lattice  imperfections  will 
again  be  neglected. 

When  one  considers  incomplete  donor  ionization  and  the 
presence  of  free  charge  carriers  within  the  depletion 
region,  the  total  space  charge  density  can  no  longer  be 
considered  a  constant.  Instead,  one  must  return  to  the 
general  expression  for  space  charge  density  given  by 

P(x)  =  e[N^(x)-n(x)+p(x)]  (4.1) 

in  which  N^(x)  is  the  density  of  ionized  donor  sites  and 
n(x)  and  p(x)  are  the  free  electron  and  hole  densities, 
respectively. 

Equation  4.1  can  be  immediately  simplified  by  assuming 
that  the  positive  space  charge  contributed  by  the  minority 
carriers  is  negligible.  This  assumption  is  valid,  since 

p(x)  «<  exp[-(Efs-Evd;/KT]  (4.2) 

and  for  N-Type  semiconductors  Efs"Evd>>KT* 

The  electron  concentration  at  the  donor  level  (nd). 
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which  will  be  used  to  calculate  N^,  may  be  obtained  by 
forming  the  product  of  the  number  of  available  energy  levels 
and  the  probability  of  occupancy.  The  donor  level  (Ed) 
actually  consists  of  narrowly  spaced  levels 5  however, 
since  the  energy  is  almost  single  valued,  the  electron 
concentration  at  the  donor  level  can  be  expressed  as 


(4.3) 


|exp[(Ed-Efs)/KT]  +  1 

where  the  term  in  the  denominator  results  from  multiplica¬ 
tion  by  the  Fermi  probability  factor,*  Since  the  electron 
concentration  at  the  donor  level  is  also  the  density  of 
unionized  impurity  atoms,  the  density  of  free  electrons 
contributed  by  ionized  impurities  must  be 


n 


c 


(4.4) 


or  by  using  Eq.  4.3, 


*The  conventional  form  of  the  Fermi  probability  factor 

for  electrons  is  f^E^)  =  [exp  [  En-EfS]  +  l]  S  however, 

\  K*  / 

the  Fermi  probability  factor  for  electrons  of  the  donor 
level  has  an  additional  factor  of  £  to  account  for  spin  degen¬ 
eracy.  For  a  rigorous  proof  of  this  form  see  S.  Wang,  Solid 
State  Electronics  (New  York:  McGraw-Hill  Book  Company,  1966), 
pp.  143-5. 
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1 

i 

I 


r* 

f 

1 


nc  ■  Nd  1  - 


*exp[(Ea-Efs)/KT]  +  1 
which  may  be  written  as 


(4.5) 


N 


nc 


2exp[-(Ed-Sfs)/KT]  +  1 


(4.6) 


Equation  4.6  may  be  expressed  in  a  more  useful  form  by  the 
following  manipulations: 

(1)  The  ionization  energy  (Ei)  is  defined  as 


Ei  =  Eco  -  Ed 


(4.?) 


(2)  By  adding  and  subtracting  Eco  in  the  numerator 
term  of  the  exponential  in  Eq.  4.6,  it  may  be  written  as 


-<Sd-Efs>  -  <Ec0-Ed>  +  Bfs  - 


CO 


(4.8) 


(3)  The  ionization  energy  can  be  incorporated  into 
Eq.  4.8  yielding 


-( E , -E_  )  =  E.  +  E_  -  E 
d  fs  i  fs  co 


so  that  an  equivalent  expression  for  Eq.  4.6  is 


(4.9) 


N, 


nc  = 


1  +  2exp(E./KT)exp[(Efs-Eco)/KT] 


(4.10) 


So 


expressed  as 


exp{E1/KT)  = 


(VNd)(1"V) 

2y2 


(4.15) 


Next,  one  can  use  the  Boltzmann  relation*  to  express 
the  actual  free  electron  density  as  a  function  of  the 
quasi-Fermi  level  for  electrons  within  the  depletion  region 
as  follows: 


n(x)  =  ncexp[-eg*(x)/KT]  (4.16) 

This  equation  represents  the  rearrangement  of  free  carriers 
in  the  conduction  band  which  occurs  to  establish  equilibrium. 
Finally,  an  expression  for  the  free  electron  density  of 
the  conduction  band  valid  in  the  depletion  region  must  be 
found.  Clearly,  this  expression  will  be  a  function  of 
distance  (x),  since  the  potential  of  an  electron  varies 
with  distance  as  shown  in  Fig.  4.1.  Eq.'s  4.6  and  4.15 
can  be  manipulated  to  show  this  dependence  and  give  the 
desired  result  by  making  the  following  substitutions: 


*The  Boltzmann  relation  states  that  the  free  carrier 
density  varies  exponentially  with  carrier  electrochemical 
potential. 


(1)  First,  one  can  write 


-(Ed-Efs)  =  -Ed  +  Efs 


(4, 


(2)  By  using  Eq.  4.7,  evaluated  in  the  depletion 
region,  the  a.bove  expression  may  be  rewritten  as 


-(Ed-Efs)  =  E 


i  +  Efs 


-  E 


cd 


(4. 


(3)  Adding  and  subtracting  Eco  yields 


)  =  E,  +  (E~  -E  )  -  (E  ,-E  )  (4. 

v  d  fs7  i  '  fs  co 7  v  cd  co7  v 


(4)  Noting  that  (£C(i-Eco)  is  by  definition  ep*(x) 
and  making  the  above  substitution  into  Eq.  k.6  yields 


f 


n  ,(x)  = 

1  +  2e:cp(Ei/KT)exp[(E,,s-Eco)/KT]exp[-ep  »/KT] 

(5)  Now  if  Eq.’s  4.11  and  4.13  are  substituted  into 
Eq.  4.20,  the  desired  result  is  obtained,  i.e., 


N. 


ncd(x)  = 


1  +  ( 1/y~1 )exp[-eP 1 (x)/KT] 


( 0=x=W ) 


(4.21) 


The  above  expression  represents  the  free  electron  density 
due  to  ionized  impurity  atoms  as  a  function  of  distance, 
valid  for  the  depletion  region  before  any  rearrangement 


X 


I 

r 

r 

r 


i 

I 


62 

occurs  due  to  electrochemical  potential.  Furthermore, 
since  each  ionized  donor  site  will  have  contributed  one 
free  electron 


(4.22) 


which  yields  the  needed  result  for  the  space  charge  density 
expression  (Eq.  4,1). 

Using  Eq.*s  4.16,  4,21,  and  4.22,  the  space  charge 
density  within  the  depletion  region  becomes 


r 


P(x)  = 


1  +  (l/Y-l)exp[-e0'(x)/KT] 
Yexp[-epf(x)/KT] 


(4.23) 


This  expression  represents  the  space  charge  density  assuming 
Incomplete  donor  ionization  and  the  presence  of  free  charge 
carriers  within  the  depletion  region.  One  should  note  that 
if  B*(x)  is  large  enough  the  term  corresponding  to  the 
free  carrier  concentration  may  be  neglected.  This  would 
correspond  to  the  Schottky  approximation  that  the  depletion 
region  is  uniformly  depleted  of  electrons.  This  appears  to 
be  a  good  approximation  for  the  region  near  x=0,  but 
significant  error  is  Introduced  in  the  region  near  x=W. 
Furthermore,  if  complete  donor  ionization  is  assumed,  i.e„, 
y=l ,  then  N^(x)  reduces  to  and  yields  the  "abrupt"  model 


which  was  the  assumed  form  for  the  Schottky  expression  for 
space  charge  density, 

Poisson* s  equation  will  now  be  used  to  relate  the  space 
charge  density  [p(x)]  to  the  negative  potential  of  an 
electron  [£*(x)3,  just  as  was  done  in  the  previous  chapter. 
However,  the  relation  will  be  considerably  more  complicated 
sJnce  both  terms  contributing  to  the  space  charge  density 
are  functions  of  distance.  Poisson* s  equation,  assuming 
a  one-dimensional  application,  becomes 

d2g  *  (x)  _  ffd £ _ 1 _ 

dx3  Ss fl  +  (l/y-l)exp[-e3* (x)/KT] 

v  (4.24) 

\ 

Yexp[-e3'(x)/KT] . 

This  expression  may  be  simplified  by  an  appropriate  change 
of  variables,  namely 

y(x)  =  ep*(x)/KT  (4.25) 


and 


u(x)  = 


e3N 


«_KT6^ 


x  =  K^x 


(4.26) 


which  upon  substitution  into  Eq.  4.24  yields. 


♦ 

i 

J 


which  is  a  second  order,  nonlinear  differential  equation. 
The  appropriate  boundary  conditions  are 


(1)  lln4|IM=0 

>:->»  dx 


(2)  0»(O)  =  Vr  +  Vd 


(4.28) 


(4.29) 


which  can  be  easily  justified  by  Fig.  4.1.  These  boundary 
conditions  must  now  be  transformed  to  the  new  variables. 
They  are 


(1)  dy/du  =  0  for  y  =  0 


(2)  y(0)  =  yQ  = 


•(vv 


(4.30) 


(4.31) 


It  is  worth  noting  that  yQ  is  a  function  of  the  applied 
bias  (Vr).  Since  Vd  is  fixed  for  a  particular  junction, 
yQ  will  in  general  be  specified  by  the  applied  bias  plus 
an  additive  constant. 

The  result  of  the  first  integration  of  Eq.  4.27  may  be 
written  as  follows; 


dy/du  =  -fz  z? 


(4.32) 


where 


z  =  ln[Yey+l-v]  +  Ye”y  -  Y 


(4.33) 
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This  expression  is  verified  in  Appendix  E. 

In  order  to  complete  the  solution,  one  must  perform  a 
numerical  integration  of  Eq.  4,32,  since  the  solution  in 
closed  form  is  not  apparent,  if  one  exists.  The  numerical 
technique  to  be  used  here  is  the  Runge-Kutta  method  of 
order  four  for  finding  a  point  by  point  solution;  utilizing 
a  digital  computer  to  perform  the  actual  calculations,  A 
discussion  of  the  Runge-Kutta  method  is  presented  in 
Appendix  F  along  with  the  computer  program  used  to  Integrate 
Eq.  4,33  and  plot  the  results.  A  typical  computer  run  is 
contained  in  Appendix  G, 

The  computer  solutions  are  shown  in  Fig.  4,2  and 
Fig,  4,3.  In  Fig.  4,2a  the  normalized  electron  potential 
y  is  plotted  against  the  normalized  distance  u.  It  has  been 
assumed  here  that  the  normalized  barrier  height  is  yo=60» 
This  would  correspond  to  a  combined  equilibrium  diffusion 
potential  and  reverse  bias  of  approximately  1.56  volts, 
i.e.,  Vd+Vr~1.56  volts.  The  degree  of  ionization  was 
varied  from  0.05  to  1,0  to  reflect  its  influence  on  the 
normalized  electron  potential  and  each  curve  corresponding 
to  a  particular  y  is  shown.  The  curves  would  seem  to 
indicate  that  the  effect  of  y  is  quite  small,  since  the 
general  shape  of  the  curves  is  maintained  as  y  is  varied. 

In  Fig.  4,2b  the  normalized  space  charge  density  (normalized 
to  unity  at  the  contact)  is  plotted  against  the  normalized 
distance  u.  For  comparison  purposes  the  equivalent  Schottky 


FIGURE  4.2a;  Computer  Solutions  for  the  Normalized  Electron 
Potential  versus  the  Normalized  Distance  into 


the  Depletion  Region.  Solution  is  for  the 
Initial  Condition  of  y  =60. 


Normalized  p 
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representation  for  space  charge  density  ( normalized  to  unity 
at  the  contact)  is  also  plotted#  The  abrupt  representation 
for  space  charge  density  remains  constant  to  a  point  wh ere 

* 

uc=W  and  V  in  terns  of  the  new  variables  y  and  u  is  W=(2yo)“. 
It  now  becomes  apparent  why  the  Schottky  model  gives  useful 
and  fairly  accurate  results#  First,  it  can  be  noted  that 
the  Schottky  and  sore  exact’  nodel  agree  to  a  point  near 
u=W={2y0)^.  After  this  point  however,  the  normalized  space 
charge  density  of  the  more  exact  nodel  decays  over  a  finite 
distance  due  to  the  increasing  compensation  offered  by  the 
increasing  concentration  of  electrons  as  x-s**#  Furthermore, 
V  (the  value  of  u  at  which  y=0)  as  predicted  by  the  sore 
exact  model  is  greater  than  that  predicted  by  the  Schottky 
model.  However,  If  one  compares  the  areas  of  disagreement 
between  the  two  models,  the  additional  space  charge  density 
predicted  by  the  more  exact  model  for  ¥>(2y0)5  has  a 
compensating  effect  on  the  reduced  s'ace  charge  density  for 
wc(2y0)^»  Thus,  the  Schottky  model  becomes  a  fairly  accurate 
model  for  an  effective  space  charge  density  and  depletion 
width* 

Similar  plots  are  shown  In  Fig#  4,3.  Here  the  value  of 
yQ  was  chosen  as  20,  to  correspond  to  a  combined  diffusion 
potential  and  reverse  bias  of  approximately  0#52  volts# 

The  Schottky  model  predicts  a  voltage  dependent 
capacitance  as  given  by  So#  3*28#  This  junction  capacitance 
arises  from  the  fact  that  the  depletion  region  acts  as  a 


Normalized  p 


Density-  versus  the  Normalized  Distance  into  the 
Depletion  Region.  Solution  is  for  the  Initial 
Condition  of  y  =20. 

v  A 
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parallel  plate  capacitor  since  for  a  small  voltage  change 
a  net  charge  is  either  added  or  substracted  from  the  region 
near  x=tf.  By  Gauss Js  law  the  total  space  charge  per  unit 
area  (QM)  contained  within  the  depletion  region  is  related 
to  the  electric  field  at  the  contact  (E0)  by 


.C*  XT 

s  o 


(4.34) 


The  electric  field  at  the  contact  is  in  turn  related  to  the 
slope  of  the  potential  function  evaluated  at  the  contact,  or 


dS  * (x) 
dx  x=0 


(4.35) 


where 


ae«(x) 

ax 


<KT/e)^j  (4y/du)|u=( 


(4.36) 


in  terns  of  the  new  variables  y  and  u.  Thus,  by  substitu¬ 
tion  Eg.  4.34  becomes 


%,  =  -  (esKTiid)^‘5y/dul  |  u=0 


(4.37) 


in  which 
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(dy/du)  u=Q  =  -/2*|ln[Yeyo  +  1  -  y]  + 


o 


(4.38) 


by  using  Ea.  4.32.  Finally,  the  junction  capacitance  can  be 
calculated  as  follows: 


r  ' 

dQ. 

w  _  v 

- I - — -  -  Ye'y° 

1  +  (l/v-l)e  y° 

dV  ~  k2 

+  Ye“yo  -  y  j 

l  J 

where 


(4.39) 


(4.40) 


yo  = 


e(V,+V  ) 
d  r 


(4.41) 


The  calculations  made  in  this  chapter  thus  far  have 
neglected  image  force  effects*  If  one  were  to  include  image 
effects,  Poisson’s  equation  (Eq.  4.24)  would  have  to  be 
modified  by  replacing  p*(r)  by  £(x)  where 


B(r)  =  g’(x)  +  £"(x) 


(4.42) 


?3 


Since  the  discussion  of  Chapter  II  indicates  that  a  single 
expression  for  g"(x),  valid  for  the  entire  depletion  region, 
is  not  possible  because  of  model  dependence  near  the  metal 
surface,  then  Eq.  4.42,  at  best,  would  have  to  be  approximated 
by  assuming  two  expressions,  one  valid  for  x>xm  and  the  other 
valid  as  x  — •>  0.  In  addition,  a  simple  transformation  of 
variables  could  no  longer  be  used  to  make  the  solutions  of 
Poisson* s  equation  Independent  of  the  particular  materials 
used  to  form  the  contact,  and  no  general  results  could  be 
obtained.  In  short,  the  addition  of  image  effects  would 
hopelessly  complicate  Eq.  4.24, 

Fortunately,  there  is  an  alternate  approach  which  can 
give  a  reasonably  accurate  approximation  to  the  maximum 
correction  of  barrier  height  needed  because  of  image  effects. 
It  has  already  been  established  that  the  Schottky  model  Is 
quite  accurate  in  the  region  near  the  junction  and  that  the 
influence  of  image  effects  predominates  in  this  region. 
Furthermore,  if  one  assumes  that  occurs  In  this  region 
and  that  xffl  is  large  enough  to  permit  (x)  to  be  approximated 
by 


isfrep 


(4.43) 


then  one  can  use  the  results  obtained  in  Chapter  III  to 
estimate  the  maximum  correction  needed  for  image  effects  by 
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E*.  In  order  to  show  the  influence  of  image  effects  on  the 
b 

result  obtained  for  the  normalized  electron  potential  as 


shown  in  Fig*  4,2a.  E£  will  be  normalized  to  y£*  This  can 
be  done  by  using  the  approximate  expression  for  E£  based  on 
the  Schottky  model,  or 


£  * 


where  EQ  is  given  by 


E 


o 


i 


(4.44) 


(4.45) 


In  terms  of  the  normalizing  variable  y,  0*(O)  may  be 
expressed  as 


0*(O) 


(4.46) 


Combining  Eq.*s  4.44,  4.45,  and  4.46,  the  normalized  correc¬ 
tion  for  image  effects  may  be  written  as 


r" 
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(e/KT) 
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L8rr2(e»e°)3€  . 

s  s  s 


■* 

■_<»o 


i 

)4 


(4.47) 


which  represents  a  "worst  case"  estimate  to  the  modification 
imposed  on  the  normalized  electron  potential  near  the  contact 
due  to  image  effects.  The  reduction  in  the  barrier  height 
yQ  is  usually  quite  small  under  normal  applications  of  bias 
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voltage  and  for  the  purpose  of  Illustrating  this  reduction 
the  following  situation  will  be  assumed: 

The  contact  is  a  metal- (N-Type) -Si  contact  with 
Vd+Vr=l #56  volts,  y=1*°»  T=300  °K,  Nd=1016  cm“3,  and 
e*rtes=12eo  farad/m.*  Under  the  assumption  that  Vd+V'r=1.56 
volts  Fig,  4,2a  is  applicable  and  one  may  proceed  to 
calculate  the  reduction  factor  y£  as  follows: 
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yb  <<  yo  =  60 


The  result  of  this  calculation  is  shown  in  Fig#  4,4, 

The  quantity  y£  seems  quite  negligible  in  comparison  to 

yQ=60  and  would  be  even  smaller  for  a  lower  doping  concen- 

16  -3 

tratJon,  i.e.,  Nd<10  cm  •  Although  this  seems  sufficient 
Justification  for  neglecting  image  effects,  one  must  be 


♦For  silicon  €g«l  since  lattice  polarization  is  purely 
21 


electronic 
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cautious  in  doing  so.  Since  y£  is  a  function  of  the  reverse 
bias  voltage  image  effects  can  be  expected  to  become  more 
important  at  higher  voltages.  In  addition,  it  has  been 
shown  in  Section  2.6  that  a  small  change  in  the  barrier 
height  would  significantly  alter  the  transport  properties 
of  the  barrier  due  to  quantum-mechanical  effects.  As  it 
turns  out,  image  effects  can  usually  be  neglected  except 
under  reasonably  large  reverse  biases,  when  the  reduction 
in  the  barrier  potential  accounts  for  the  voltage  dependence 
of  the  reverse  saturation  current.  Image  effects  will  be 
discussed  again  in  relation  to  the  I-V  characteristic  of 
the  contact  in  Chapter  VI. 
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CHAPTER  V 

METHODS  0?  DSTSRRIiSBiG  THE  EQUILIBRIUM  BARRIER 
HEIGHTS  0?  KS2AL-  {  N-TYFS  )  -SSKICCJsBDCTOR 
CONTACTS:  A  CCHPARISC^  C?  THE  MODELS 
CF  CHAPTER  III  AND  CHAFER  IV 

The  most  commonly  used  experimental  methods  to  determine 
the  energy  barrier  heights  and  evt  are  by  the  photo- 
re  spoase  xethod  and  by  an  extrapolation  of  data  obtained 
from  a  capacitance  versus  voltage  plot.  In  this  chapter 
both  methods  will  be  discussed;  however ,  major  attention 
will  be  focused  on  the  capacitance  method  since  the  results 
obtained  are  dependent  on  the  particular  model  used  for  the 
metal- semiconductor  contact*  Furthermore,  since  the  photo- 
response  method  provides  a  direct  measure  of  the  barrier 
height,  a  comparison  of  the  results  predicted  by  the  two 
methods  will  provide  a  means  of  comparing  the  accuracy  of 
the  contact  models  discussed  in  Chapters  III  and  IV. 

The  photoresponse  method  uses  a  monochromatic  light 
source  to  induce  a  photocurrent  in  the  device  under  investi¬ 
gation.  The  source  may  be  incident  upon  either  the  semi¬ 
conductor  or  metal  side  of  the  contact.  To  provide  a  direct 
measurement  of  eV^  the  source  must  be  incident  upon  the 
semiconductor  and  to  measure  the  source  must  be  incident 
upon  the  metal;  in  either  case,  the  material  upon  which  the 
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light  source  is  incident  must  be  as  thin  as  possible  to 
enable  sufficient  light  for  electron  excitation  in  the 
immediate  vicinity  of  the  contact.  Generally,  the  light 
source  is  incident  cn  the  re  tel  since  the  mc-st  common  method 
of  preparing  ratal- seal  conductor  contacts  is  by  vacuus 
evaporation  of  the  setel  onto  a  doped,  single  crystal 
substrate*  Sals  method  allox ^  very  accurate  control  of  the 
metal  thickness  and  evaporation  of  a  r&tel  has  the  addi¬ 
tional  advantages  of  lexer  evaporation  temperature s  while 
avoiding  the  problems  of  doping  and  crystal  structure. 

Assuming  that  the  light  source  is  incident  on  the  metal. 

Fig*  5*1  def  ines  the  relevant  parameters  involved  in 
measuring  E-0  by  the  photorespanse  method* 

Mhen  light  is  directed  onto  the  metal  surface  the  free 
electrons  in  the  metal  are  excited  and  if  sufficient  energy 
is  acquired,  they  will  overcome  the  energy  barrier  (E^) 
and  produce  a  current  in  the  semiconductor*  if  the  photo¬ 
electric  current  per  incident  photon  is  defined  as  the 
photoresponse  (H),  then  the  spectral  distribution  of  E 
will  reveal  the  barrier  height  (Efe)  as  seen  by  an  electron 
leaving  the  metal*  A  typical  spectral  response  is  shown 
in  Fig*  5*2,  Figure  5*2  also  shows  that  there  are  two 
distinct  regions  of  the  photoexcitation  process: 

(1)  The  photoemission  of  excited  electrons  in  the 
metal  over  the  barrier  E^. 

(2)  The  band  to  band  excitation  of  electron-hole  pairs 
’ithin  a  diffusion  length  of  the  semiconductor  depletion  region* 


FICUSS  5,1;  Energy  Diagram  of  a  Eetal-Seslccnducfcor  acetifying 
Contact  Shoving  the  Parameters  Involved  In  the 
Measurement  of  E.  by  the  Photoexcitation  Hethod. 
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■The  Foxier  theory2^  of  the  photoelectric  effect  predicts 
that  the  dependence  of  the  photocurrent  { J^)  on  photon 
energy,  for  photon  energies  exceeding  S^+3£T*  r»y  be 
expresceu  as5* 

J  ac  (fcf  -  hf„}a  (5*1) 

o  c 


where 


(5.2) 


as  shown  in  Fig,  *>,1 »  Thus,  if  one  plots  2^  (which  is 
defined  to  be  proportional  to  Jp)  versus  the  incident 
photon  energy  (hf),*  a  straight  line  would  result,  which  in 
turn  yields  an  extrapolated  value  for  E^,  i.e.,  the  intercept 
value  on  the  hf  axis*  A  plot  of  2^  versus  hf  for  an 
Au-CdSe  contact  is  shown  in  Fig,  5,3,  Figure  5, J  also 
illustrates  quite  clearly  why  it  is  necessary  to  have  a  very 
thin  metal  contact,  since  an  accurate  extrapolation  for 
becocjes  difficult  for  metal  widths  exceeding  900  2, 

The  capacitance  versus  voltage  technique  for  determining 


*The  values  of  R  are  taken  from  the  longer  wavelengths 
of  f  to  insure  data  corresponding  to  the  process  of  photo- 
emission  from  the  metal  and  not  band  to  band  excitation. 


I 

i 


|S~  (arbitrary  units) 


Photoresponse  of  a  Au-CdSe  Contact  versus  Fhoton 
Energy?^  The  Curve  Labeled  A  if  for  a  Ketal  Thick¬ 
ness  of  300  R  and  the  Curve  Labeled  B  is  for  a  Metal 
Thickness  of  900  R. 
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the  equilibria  diffusion  potential  (V^) ,  assuring  an 
idealized  contact,  can  also  provide  an  indirect  seasure- 
cent  of  S^j,  since  Eb  is  related  to  as  follows: 


(5.3) 


which  is  evident  frcn  Fig*  2.5  with  V&=0.  The  quantity 
vTd  is  experimentally  determined  by  plotting  the  square  of 
the  reciprocal  of  the  Junction  capacitance  per  unit  area 
(C  j)  versus  the  reverse  bias  voltage  (Vr).  A  plot  of  this 
fora  will  allow  one  to  extrapolate  a  value  of  (the 
intercept  value  on  the  Vr  axis)  and  once  V.  is  known,  E^ 
rsay  be  calculated  by  using  Eq®  5*3*  The  najor  difficulty 
encountered  when  using  this  method  is  that  the  theoretical 
expression  for  the  Junction  capacitance  (Eq. 1  s  3*28  and  4.39) 
Is  "model  dependent”  and  the  relation  of  the  experinentally 
deterained  intercept  voltage  to  Vd  will  reflect  this 
dependence.  This  will  be core  apparent  in  the  discussion 
which  follows. 

First,  consider  the  capacitance  expression  based  cn 
the  Schottky  aodel  (Eq.  3*28).  If  this  expression  is  squared 
and  the  reciprocal  taken  of  the  result,  one  obtains 


2 (V.+V  ) 
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Next,  consider  the  capacitance  expression  derived  on  the 


basis  of  the  more  exact  model  of  Chapter  IV,  Squaring 
£q,  4.39  and  taking  the  reciprocal  yields 

(5.10) 


'■f  =  E* 


j  jlnLYe^°+l-Y]  +  ye“^°  -  y 


3  f - 1 - =  -  Ye-y°T 

l\l  +  (l/y-l)e  y°  }) 


in  which 


b 


2KT 


(5.11) 


This  expression  Is  considerably  more  complicated  than  the 
equivalent  expression  (Eq.  5»*0  based  on  the  Schottky  model. 
However,  for  sufficient  values  of  reverse  bias  yQ»l, 

Eq.  5.10  may  be  approximated  by 


’i 


»  k^‘’^ln[Yey°J*l-Y]  - 


(5.12) 


Expanding  the  logrithmic  term  by  the  power  series 
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a  further  approximation  may  be  made  for  yQ»l , 
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yielding* 


c“2  »  k~2[yQ  +  ln(y)  -  Y]  (5.14) 

one  may  now  recognize  the  equation  of  a  straight  line  if 
(k^/Cj)2  is  plotted  against  yQ.  The  approximations  leading 
to  Eq.  5.14  can  be  easily  justified  by  observing  a  plot 
of  Eq.  5.10  for  arbitrary  values  of  yQ  and  y.  Such  a  plot 
is  shown  in  Fig.  5.5.  For  values  of  yQ>4  (approximately), 

Eq.  5*10  becomes  a  straight  line  and  an  extrapolation  of 
th^se  lines  for  each  y  may  be  used  to  determine  the  intercept 
value  of  VQ.  This  intercept  value  of  VQ  is  of  primary 
importance  since  it  will  result  in  an  expression  for  Vd. 

The  quantity  Vd  can  be  obtained  by  noting  that  when 
(k-^/Cj)  =0,  the  straight  line  approximation  to  Eq.  5.10 
can  be  written  as 


27 

*Dewald  '  has  obtained  the  same  limiting  result  for 
C j2;  however,  his  original  expression  for  C "j2  is  not  the 
same  as  Eq.  5*10  although  both  expressions  produce  the  same 
limiting  result  as  reverse  bias  becomes  large.  Further¬ 
more,  his  derivation  of  a  voltage  dependent  capacitance 
follows  from  an  investigation  of  the  distributions  of 
charge  and  potential  at  a  zinc  oxide-electrolyte  inter¬ 
face. 


Bias  Voltage  (Eq.  5»t0)„ 


(5.15) 


:»o  T  ln(y)  -  Y  =  0 
or 

e(V-7/) 

— »  y  -  ta(Y)  (5.16) 

Rearranging  yields 

Vd  ~  Vo  +  (5T/e>CY-ln(Y)]  (5.1?) 

It  should  be  noted  that  the  above  expression  for  the 
equilibrium  diffusion  potential  (?d)  has  neglected  the 
influence  of  Image  effects.  If  image  effects  were  included, 
then  yQ  in  Sq.  5* 15  would  be  replaced  by  the  effective 
barrier  height  near  the  contact,  I.e-,  yo-y£.  However, 
since  y0»y£  one  would  seem  justified  in  neglecting  y£. 

For  this  particular  application  the  error  introduced  by 
neglecting  y£  is  quite  small  and  It  will  become  even  smaller 
as  the  bias  voltage  is  increased.  This  may  not  be  the  case 
when  considering  the  carrier  transport  properties  of  the 
contact  since  the  magnitude  of  y£  increases  as  the  bias 
voltage  is  increased  and  quantum-mechanical  effects  are 
critically  dependent  on  the  barrier  height  and  the  electric 
field  strength  at  the  contact,  which  also  increases  with 
bias.  Thus,  one  can  safely  neglect  image  effects  for  this 
particular  application  and  Eq.  5»17  should  give  a  reasonably 


accurate  relation  between  the  intercept  voltage  and  the 
equilibrium  diffusion  potential;  however,  when  one  con¬ 
siders  the  current- voltage  properties  of  the  contact  inage 
effects  nay  become  fairly  inportant  at  higher  reverse  bias. 

The  slope  of  the  straight  line  approximation  for 
Eq.  5»10  reduces  to  the  sane  slope  as  obtained  using  the 
Schottky  expression  for  junction  capacitance.  This  can  be 
easily  verified  by  taking  the  derivative  of  Eq.  5»1^  with 
respect  to  voltage.  In  other  words,  the  value  for 
predicted  by  the  two  nodels  is  exactly  the  same  (Eq.  5*9) 
since  the  slopes  agree  for  yQ»l.  This  is  shown  quite 
clearly  in  Fig*  5.5. 

By  comparison  of  Eq.*s  5»?  and  5-1?  one  can  see  that 
there  will  be  a  significant  difference  between  the  diffusion 
potential  predicted  by  the  Schottky  and  more  accurate  model 
when  using  the  capacitance  versus  voltage  method  to 
determine  VQo  The  exact  disagreement  between  the  two  models 
will  be  the  difference  between  the  intercept  values  (V0) 
as  shown  in  Fig.  5.5.  and  will  be  Influenced  by  the  degree 
of  ionization.  Thus,  one  can  expect  the  value  of 
predicted  by  the  Schottky  model  tc  be  in  error  by  at  least 
KT/e  (from  Eq.  5»1 7  with  y=1)  and  this  error  would  Increase 
as  the  degree  of  ionization  decreases.  Furthermore,  the 
value  of  the  doping  concentration  (N^)  obtained  by  the 
capacitance  method  will  be  independent  of  the  model  used, 
since  the  slope  predicted  by  the  Schottky  model  and  the 
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straight  line  approximation  to  Eq.  5*10  are  the  sane  for 
sufficient  values  of  reverse  bias. 

The  conclusions  drawn  thus  far  are  for  an  idealized 
contact  and  one  will  find  that  an  application  of  the 
capacitance  method  to  determine  Vg  for  a  "real"  contact 
will  involve  many  complications.  Goodman  has  made  a 
comprehensive  study  of  the  complications  Involved  in  the 
measurement  of  barrier  heights  by  the  capacitance  method  for 
"real"  metal- semiconductor  contacts,  and  outlines  a 
procedure  to  minimize  many  of  these  complications.  He 
also  examines  the  model  dependence  of  the  capacitance 
expression  and  concludes  that  "the  true  value  of  the  barrier 
height  is  greater  than  the  ‘intercept  value*  by  KT/e." 

However,  there  can  be  significant  error  in  this  conclusion 
for  Incomplete  ionization,  i.e.,  y<1*  as  can  seen  from 
Eq.  5.17. 

Having  determined  the  diffusion  potential  by  experimental¬ 
ly  determining  VQf  one  can  calculate  the  barrier  height  as 
seen  by  electrons  of  the  metal  (Eb)  by  using  Eq.  5*3.  In 
order  to  use  Eq.  5*3  the  position  of  the  Fermi  level  must 
be  known  for  the  bulk  semiconductor  and  this  may  be  calculated 
by  using  the  following  relation:2^ 


i( 
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KTln 


(5.18) 
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Equation  5*18  is  an  exact  expression  and  all  of  the 
quantities  are  either  known  or  material  constants,  with  the 
possible  exception  of  which  can  be  determined  by 
Ea,  5*9-  Thus,  using  Ea.*s  5» 3.  5»17  and  5<>l8  (V^  by  the 
more  exact  model)  or  5«7  (V,,  by  the  Schottky  model)  &  value 
of  Sjj  can  be  determined  and  a  comparison  can  be  made 
between  as  measured  by  the  photoresponse  method  and 
indirectly  measured  by  the  capacitance  method.  This 
comparison  would  be  an  indication  of  the  degree  of  accuracy 
associated  with  the  Schottky  and  more  exact  models; 
however,  any  comparison  should  be  made  with  reservation 

OO 

since  Goodman  has  indicated  that  there  may  be  many 
experimental  errors  associated  with  experimentally  determin¬ 
ing  VQ. 
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CHAPTER  VI 

CURRENT-VOLTAGE  CHARACTERISTICS  OF  A  KETAL- 
S3NIC0NDUCT0R  RECTIFYING  CONTACT 

6.1  Current-Voltage  Equation. 

Before  undertaking  the  development  of  the  current- 
voltage  equation  of  a  metal-semiconductor  contact  it  would 
be  helpful  to  understand  more  clearly  the  basic  mechanism 
by  which  rectification  takes  place.  In  the  ab fence  of  an 
external  voltage  (Fig.  2.4d)  the  electrons  of  the  metal  are 
in  dynamic  equilibrium  with  the  conduction  electrons  of  the 
semiconductor.  Under  this  condition  the  rates  at  which 
carriers  transverse  the  barrier  are  equal  from  either 
direction  and  the  probability  of  such  a  crossing  taking 
place  will  depend  on  the  number  of  electrons  having  thermal 
energies  greater  than  Eb  and  moving  in  the  right  direction. 
Furthermore,  this  rate  of  carrier  transfer  will  be  an 
exponential  function  of  the  barrier  height  and  would  thus 
decrease  rapidly  as  the  barrier  height  is  increased.  Since 
the  probability  for  electron  transfer  is  the  same  for  either 
direction,  one  would  expect  no  net  electron  current  to  flow 
across  the  junction. 

Next,  consider  the  situation  in  which  an  external  bias 
is  introduced  in  such  a  way  as  to  make  the  semiconductor 
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positive  with  respect  to  the  metal  (Fie;.  2.5b).  This  would 


correspond  to  a  reverse  bias,  since  the  effective  barrier 
height  to  electron  flow  from  the  semiconductor  into  the  metal 
has  been  increased.  On  a  probability  basis,  the  probability 
of  electron  transfer  from  the  metal  is  still  governed  by 
exp(-Eb/KT)  and  must  therefore  remain  ’unchanged  with  respect 
to  the  equilibrium  condition.  On  the  other  hand,  the 
probability  of  electron  transfer  from  the  conduction  band  of 
the  semiconductor  into  the  metal  will  be  proportional  to 
exp[-e(V^+Vr)/KT3  and  must  therefore  be  greatly  reduced  even 
for  a  small  amount  of  reverse  bias  (Vr).  The  equilibrium 
condition  of  equal  electron  transfer  can  no  longer  be  main¬ 
tained  and  a  small  net  current  will  cross  the  Junction.  This 
will  be  called  the  reverse  saturation  current  (IQ)  and 
represents  an  electron  leakage  current  from  the  metal  into 
the  semiconductor  under  reverse  bias.  Here  the  term 
"leakage"  is  used  to  emphasize  the  fact  that  it  is  an 
unwanted  current,  since  Ideally  I  should  be  zero  for  perfect 
rectification. 

Finally,  consider  the  situation  in  which  the  polarity 
of  the  bias  voltage  is  reversed,  i.e.,  the  metal  is  now  made 
positive  with  respect  to  the  semiconductor.  This  will 
correspond  to  a  forward  biased  condition  and  is  shown  in 
Fig.  2.5a.  Again  equilibrium  conditions  cannot  be  maintained, 
since  the  probability  of  electron  transfer  from  the  conduc¬ 
tion  band  of  the  semiconductor  to  the  metal  is  greatly 
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increased*  i*e.t  it  must  now  be  proportional  to 
eip[-e(V^-V^)/KT].  Furthermore,  the  rate  at  which  electrons 
flow  in  the  opposite  direction  remains  unchanged  since  this 
transfer  is  still  proportional  to  expC-E^/KT)  and 
remains  constant  (neglecting  image  effects).  The  result  of 
this  unbalance  is  a  large  net  electron  current  flowing  into 
the  metal  under  the  influence  of  a  forward  bias. 

Thus,  the  character  of  rectification  is  exhibited 
through  the  unbalanced  flow  of  electron  current  under  forward 
and  reverse  bias  conditions.  In  addition,  one  should  note 
that  for  a  metal- (N-Type)- semiconductor  contact  positive 
values  of  V_  correspond  to  a  reverse  bias  and  negative 

a 

values  of  V„  to  a  forward  bias.  Furthermore,  the  current 

a 

flowing  under  forward  bias  will  be  of  much  greater  magnitude 
than  under  reverse  bias  conditions. 

In  order  to  calculate  an  expression  which  can  predict 
the  current-voltage  relationship  of  a  metal-semiconductor 
contact  it  is  necessary  to  assume  that  the  actual  number  of 
electrons  constituting  current  flow  across  the  contact  is 
only  a  small  fraction  of  the  total  free  electron  population 
of  the  semiconductor.  This  seems  to  be  a  reasonable  assump¬ 
tion  for  effective  barriers,  i.e.,  Ejj»KT,  and  allows  one 
to  assume  that  the  electron  concentration  of  the  bulk 
semiconductor  is  constant  and  Independent  of  the  current 
flowing  in  the  device.  Furthermore,  it  is  possible  to 
distinguish  two  types  of  models  for  rectification  on  the 
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basis  of  barrier  width.  The  most  general  case  is  to  consider 
the  width  of  the  barrier  within  an  energy  increment  of  KT 
near  the  top  as  being  large  compared  with  the  mean  free 
path*  (Ae)  of  an  electron  in  the  semiconductor.  This 
distance  near  the  top  of  the  barrier  is  used  because  it  is 
necessary  to  compare  the  probability  of  an  electron  being 
stopped  by  a  normal  collision  process  and  the  probability 
of  its  reflection  by  the  potential  barrier  itself.  When  the 
width  within  KT  of  the  top  of  the  barrier  is  large  compared 
to  the  mean  free  path  of  an  electron,  the  electrons  crossing 
the  barrier  can  be  expected  to  suffer  many  collisions  before 
reaching  the  other  side  and  as  a  consequence  current  flow 
across  the  barrier  must  be  by  a  diffusion  process.  When 
this  width  is  small  compared  to  Xe,  the  probability  of  many 
collisions  taking  place  is  small  and  current  flow  can  be 

30 

thought  of  as  an  emission  of  electrons  over  the  barrier. y 
The  emission  model  is  somewhat  simpler  than  the  diffusion 
model,  but  the  diffusion  model  must  be  regarded  as  the  most 
general  since  Xe  is  on  the  order  of  10  5  to  10“  cm, 
causing  the  model  to  be  applicable  to  all  but  the  thinnest 
barriers.  A  possible  example  of  a  thin  barrier  might  be  a 
point  contact  device,  in  which  case  xo  could  exceed  the 

V 

*\  is  defined  as  the  average  distance  between  two 
successive  collisions  of  an  electron  with  the  lattice  struc¬ 
ture  of  the  semiconductor. 
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width  of  the  barrier  near  the  top. 

Under  the  assumption  that  diffusion  theory  is  applicable 
and  that  the  potential  within  the  barrier  is  known  from  the 
model  developed  in  Chapter  IV,  one  may  calculate  the 
current- voltage  relationship  for  the  contact.  Quantum- 
mechanical  effects  and  image  force  will  be  neglected 
initially  so  that  the  current  flow  will  be  assumed  to  be 
strictly  by  a  diffusion  process.  Furthermore,  since  the 
current  in  the  barrier  region  depends  on  the  local  field  as 
well  as  the  local  concentration  gradient,  a  calculation  of 
the  current-voltage  relation  for  a  N-Type  material  must 
begin  from  the  general  diffusion  equation  of  the  fornr 


..rv  Tdn(x)  ,  eE(x)n(x) 
8A3nrto  +  KT 


(6.1.1) 


The  diffusion  of  carriers  against  a  retarding  electric 
field  will  occur  across  the  barrier  if  a  difference  in 
electron  density  occurs,  and  with  the  aid  of  Fig,  6,1  this 
difference  in  respective  densities  can  be  easily  shown.  The 
minimum  value  of  electron  concentration  will  ocur  at  the 
junction,  i.e.g 


n(  0}  =  Y^ciexp( -eV^/KT)  (6.1.2) 

One  should  note  that  this  density  will  be  independent  of 
the  applied  voltage  since  the  electron  concentration  in  the 


Barrier  Under  Forward  Bias  Conditions; 
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metal  and  the  barrier  height  have  constant  values 
(neglecting  image  effects).  Furthermore,  n(x)  will  increase 
as  x  increases  until  it  reaches  a  maximum  at  x=W  where 


n(W)  =  no  =  yNd 


(6.1.3) 


and  may  also  be  regarded  as  constant  since  it  was  originally 

assumed  that  the  equilibrium  concentration  of  electrons 

in  the  bulk  semiconductor  is  not  appreciably  disturbed  by 

current  flow  in  the  device.  Thus,  a  difference  in  electron 

densities  between  opposite  sides  of  the  barrier  is  apparent 

and  one  can  expect  a  diffusion  of  electrons  in  the  direction 

of  the  metal  through  the  potential  barrier.  Electron 

diffusion  constitutes  an  electron  current  (I  )  and  if  one 

n 

adopts  the  convention  that  positive  current  flows  in  the 
negative  x  direction  the  calculation  of  In  may  proceed  as 
follows: 

(1)  Equation  6.1.1  is  multiplied  by  exp[-eg’ (x)/KT] 
and  upon  substitution  of  E(x)=-dg*  (x)/dx  yields 


(x)/KT]  =  eADn 


dn(x) 

dx 


(6.1.4) 


e  .  dg 1 (x) 
KT  dx 


n(x) 


exp[-eg'(x)/KT] 


(2)  Noting  that 
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^  exp[-ef3‘(x)/KT]  = 


(6.1.5) 

-  *  exp[-eg'(x)/KT] 


the  right  side  of  Eq.  6.1.4  may  be  expressed  in  terms  of  a 
derivative,  or 


In©xp[-e3*(x)/KT]  = 


eADn'  ^|^n(x)exP[-eP  1  (x)/KT]J 


(6.1.6) 


(3)  Multiplying  both  sides  of  Eq.  6.1.6  by  dx  and 
integrating  over  the  depletion  width  yields 


J 


f  W 

I^xpC-ePMxVKT]  = 


eADn^n(x)exp[-eg' (x)/KT]J 


(6.1.7) 


(4)  The  right  side  of  Eq.  6,1,7  may  be  evaluated  by 
noting  the  following  boundary  conditions* 


(a)  n(W)exp[-eg’ (W)/KT]  =  YNdexp[-e(V(i-V;)  )/KT]  (6,1.8) 
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(fc)  n(0)expOe3'(0)/KT]  =  YNdexp[-eVd/KT]  (6.1.9) 

(5)  Finally,  I  may  be  regarded  as  a  constant  if 
recombination  is  neglected  in  the  depletion  region.  Thus, 
zaay  be  expressed  as 


*n  = 


eAQnyNdex?[ -eVd/KT][exp(ey ^/KT?  -  l] 


Jr 


(6.1.10) 


exp[ -ep ' lx) /KT]dx 


Equation  6.1.10  takes  the  form  of  the  familiar  diode 
equation  and  may  be  rewritten  in  a  shortened  form  as 
f oi lows : 


T 

~n 


(6.1.11) 


where 


I 

o 


eA D^y ’>d  J^exp  ( -e Vd/KT  Jj  • 

(Jr  [-eg«(x)/KT]j 


(6.1.12) 


Upon  closer  examination  of  the  integral  which  appears  in 
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Eq.  6*1 *1 2 ,  one  finds  that  after  the  variable  transformation 
used  in  Chapter  IV  (Sq.*s  4,25  and  4.26)  I  may  be  expressed 


(6,1.13) 


and  the  integral  could  be  evaluated  by  using  numerical 
techniques  similar  to  those  used  to  evaluate  y.  This  is 

•y 

pointed  out  because  a  numerical  evaluation  of  /e*au  would 
be  a  simple  matter  of  inserting  several  statements  into  the 
computer  program  which  has  already  been  used  to  calculate 
the  normalized  electron  potential  and  space  charge  density 
(Appendix  ?).  However,  since  each  evaluation  of  /eydu  would 
depend  on  the  Initial  value  of  y  (a  function  of  bias)  and 
the  degree  of  ionization  (Eq.  D.7),  one  must  conclude  that 
Sa.  6,1.13,  and  thus  Eq.  6.1.12  already  appear  in  their 
most  general  form. 

Thus,  If  one  neglects  quantum-me chan i cal  and  Image 
effects,  the  I-V  characteristic  may  be  stated  in  its  most 
general  form  as  follows: 


I  =  I  [exp(eV  /KT )  -  l] 
n  o  a  j 


(6.1.14) 


where  IQ  is  given  by  either  Eq.  6.1.13  (in  normalized  form) 
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or  by  Eq.  6.1. 12*  Inherent  in  Eq.  6.1.14  is  the  assumption 
that  the  magnitude  and  polarity  of  V&  is  consistent  with  the 
assumptions  on  which  the  diffusion  model  is  based. 


6.2  Influence  of  Tunneling  and  Qaantum-Kechanical 
Reflection  on  the  I-V  Characteristic. 

In  the  previous  section  quantum-mechanical  effects 
were  neglected  in  forming  the  diffusion  model  for  current 
transport  at  a  metal- semiconductor  contact.  Clearly,  the 
discussion  as  presented  In  Section  2.6  would  Indicate  that 
the  effects  of  tunneling  and  QKR  may  exert  considerable 
influence  on  the  overall  I-V  characteristic  of  the  contact. 

The  effects  of  tunneling  and  QKR  can  be  Incorporated 
into  the  diffusion  model  by  introducing  a  factor  f^  which 
represents  the  ratio  of  total  current  flow  (1^.)  predicted 
considering  quantum-mechanical  effects  to  the  diffusive 
current  flow  neglecting  these  effects  (In).  Stated 
mathematically , 


It  =  f  I 
t  q  n 


(6.2.1) 


Furthermore,  a  numerical  value  for  fq  may  be  calculated  by 
averaging  the  QKIC  over  the  complete  kinetic  energy  spectrum 
associated  with  the  electrons  incident  on  the  potential 
barrier  (corrected  for  image  effects).  If  the  kinetic 
energy  associated  with  an  Incident  electron  is  denoted  by 
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and  it  is  assumed  that  EQ  is  governed  by  a  Maxwellian 
distribution,  then  may  be  computed  as  follows: 


MTC ) exp ( -E  /KT )  dE 

6  v 


(6.2.2) 


Crowell  and  Sze'?-J  have  used  numerical  techniques  to 
compute  fa  for  Au-GaAs,  Au-Ge,  and  Au-Si  N-Type  contacts 
as  a  function  of  the  electric  field  at  the  contact  (EQ) 
for  selected  temperatures.*  A  plot  of  fQ  versus  EQ  is 
shown  in  Fig.  6.2  using  the  curves  of  the  QKTG  shown  In 
Fig.  2.9  for  an  Au-GaAs  (N-Type)  contact  for  selected 
temperatures.  For  this  particular  contact  quantum- 
mechanical  effects  may  be  seen  to  exert  considerable 
influence  on  the  total  current  flow  as  predicted  by  the 
diffusion  model  since  tunneling  becomes  excessive  for 
electric  fields  exceeding  10^  v/cm.  Similar  results  are 
found  for  the  Au-Sl  contact  over  the  same  range  of  Ec; 
however,  for  the  Au-Ge  contact  f  approaches  a  very  low 
value  at  Eo>10"'  v/cm  (fQ«0.01 ) .  This  would  indicate  that 
reflection  predominates  in  the  germanium  contact  at  relatively 
high  electric  fields. 

An  alternate  approach  to  calculating  for  a  contact 


*They  h^ve  assumed  the  approximate  form  of  E£  as 
”  eSr,l5 

• 

C" 
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barrier  would  be  to  adopt  the  approximate  expressions  given 
19 

by  Kemble  (Appendix  C)  for  the  QMTC  and  aveiage  these 
approximate  expressions  assuming  a  Maxwellian  distribution 
in  electron  energy.  However,  a  calculation  cf  this  type 
would  also  require  numerical  techniques  because  of  the 
complexity  of  the  expressions  for  the  QMTC. 

Still  another  approach  is  available  if  one  is  willing 
to  approximate  a  QMTC  plot  by  simple  analytic  expressions 
in  order  to  allow  the  average  value  of  f  to  appear  in 
closed  form.  In  order  to  illustrate  this  technique  the 
QMTC  will  be  approximated  by  a  single  straight  line  whose 
slope  has  been  adjusted  to  obtain  the  best  agreement  possible. 
A  plot  of  this  form  is  shown  in  Fig.  6.3.  This  would 
certainly  represent  a  crude  approximation,  but  Increased 
accuracy  would  be  obtained  by  using  additional  straight 
lines,  higher  order  analytic  expressions,  or  a  combination 
of  both.  Using  the  single  straight  line  approximation  the 
calculation  of  f^  would  proceed  as  follows: 

(1)  f  =  ( QMTC ) exp ( -E  /KT ) dE  (6.2.3) 

A  e  e 

(2)  In  terms  of  the  straight  line  approximation  for 
the  QMC,  f  may  be  expressed  as 


fa  ~  KT 


f  E] 

?J(0,e 


7* 

E  +  E 

(0)exp(-Ee/KT)dEe  + 


< 


-/ 

-E, 


exp( 


(S.2.*l) 


-Ee/KT ) dEe  +  ^ J ( 1 ) exp ( -Ee/KT ) dEe 


(3)  Integration  with  resoect  to  E  and  substitution  of 
the  limits  yields 


KT 


E1  +  Er 


exp(E1/KT) 


-  exp(-Er/KT) 


) 


(6.2,5) 


(4)  Since  plots  of  the  QMTC  reveal  that  E  >E. ,  a 

r  ± 

reasonable  approximation  would  be 


_  KTexp(E1/XT) 


E1  +  E; 


(6.2.6) 


Thus,  using  this  form  of  approximation  for  the  QHTC,  f 

may  be  expressed  in  terms  of  the  end-point  energies  Er 

and  E^.  One  should  also  note  that  these  end-point  energies 

are  rather  strong  functions  of  Eo,  and  thus  a  function  of 

the  applied  bias.  In  addition,  if  more  accurate  approximations 

are  used  for  the  QKTC  then  would  become  a  more  complicated 

expression  involving  more  characteristic  energies. 


6.3  Influence  of  Image  Effects  on  the  1-7  Characteristic. 


The  sample  calculation  presented  at  the  end  of  Chapter 
IV  indicated  that  the  normalized  correction  to  the  barrier 


no 


height  due  to  image  effects  (y£)  was  quite  small  in  com¬ 
parison  to  the  uncorrected,  normalized  barrier  height  at 
the  contact  (yQ).  However,  the  importance  of  image  effects 
should  not  be  overlooked  since  a  small  change  in  the  barrier 
height  has  significant  influence  on  the  QMTC.  This  can  be 
seen  from  Fig.  2.9  in  which  the  QKTC  changed  significantly 
for  small  changes  in  the  kinetic  energy  associated  with  an 
incident  electron  when  the  kinetic  energy  is  referenced  to 
the  top  of  the  barrier  corrected  for  image  effects. 
Furthermore,  y£  is  voltage  dependent  so  that  greater  influ¬ 
ence  may  be  expected  at  higher  bias  voltages.  As  it  turns 
out  y£  exerts  greater  Influence  at  higher  reverse  bias, 
but  due  to  large  currents  flowing  under  forward  bias 
(eeVf/KT»i)  the  Influence  of  y£  may  go  completely  un¬ 
noticed. 

In  order  to  include  the  influence  of  image  effects  in 
the  I-V  characteristic  predicted  by  the  diffusion  model 
the  original  assumption  that  n(0)  was  independent  of  bias 
(assuming  Efc  was  independent  of  bias,  Eq.  6.1.12)  must  be 
modified  to  include  the  Influence  of  image  effects.  This 
can  be  done  by  replacing  Vd  by  an  effective  diffusion 
potential  (VdsQ)  which  Includes  image  effects.  Stated 
mathematically. 


V,  •=  V,  -  a.  (V  ) 
dso  d  b  a 


(6.3.1) 


ill 


in  which  afc(Va)=E* (Va)/e  and  represents  the  voltage  dependent 
reduction  in  the  diffusion  potential  due  to  image  effects. 
Furthermore,  ab(Va)  may  be  written  in  the  form 


ab(V  '  %oab 


(6.3.2) 


in  which  a^Q  is  independent  of  bias  voltage  and  contains 
the  voltage  dependence.  Expressions  for  a^0  and  a£  may  be 
found  by  writing  Eq.  3*3.4  in  the  form* 


r  e3Nd  T  i 

b  a  [_8n2(efeo)2€  J  d 

S  S  5 


(6.3.3) 


so  that  ahri  may  be  expressed  as 


a.  = 

oo 


e3t!dvd 

8n3(€'e»)se 

O  O 


(6.3.4) 


*b<V  -  (1  "  W* 


(6.3.5) 


It  should  also  be  noted  that  the  inclusion  of  image 
effects  will  alter  the  integral  expression  in  Eq»  6,1.12 
since  £‘{x)  must  now  be  replaced  by  the  total  electron 


♦Recalling  that  Va=-Vr  for  reverse  bias  and  Va«Vf 
for  small  forward  bias. 
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potential  g(x)  where  g(x)  is  given  by 

g(x)  =  g'(x)  -  (x>0)  (6.3.6) 

s 

6.4  Complete  I-V  Characteristic  for  a  Metal- 
Semiconductor  Contact. 

The  complete  expression  for  the  I-V  characteristic  for 
an  idealized  metal-semiconductor  contact,  assuming  the 
diffusion  model  is  applicable,  may  be  summarized  as 
follows: 

It  =  Pexp(-eVH/KT)[exp(eVa/KT)  -  l]  (6.4.1) 


where 


eAD  vN , f  exp( -a,  a '/KT) 

T  t  _  p  d  q  bo  b  , 

o  " 


I 


(6.4.2) 


exp[~eg (x)]dx 


i 


fq  =  KT  I  (Q^‘TC)exp(-Ee/KT)dEe 


(6.4.3) 


'bo  = 


:  f%vd 


(6.4.4) 
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(1  - 


w* 


(6.4.5) 


p(x)  =  e»(x)  -  . 

s 


(6.4.6) 


and 


3f(x}  is  the  solution  to  Eq.  4.24  (6.4.7) 

6.5  Breakdown  Mechanisms  for  a  Metal-Semiconductor 
Contact. 

At  large  applications  of  reverse  bias,  junction 
breakdown  may  occur  by  anyone  of  the  following  mechanisms: 

(1)  The  barrier  becomes  so  thin  and  the  electric 
field  so  large  that  current  flow  due  to  tunneling  may  become 
excessive. 

(2)  Electrons  can  travel  to  the  conduction  band  by 
tunneling  from  deep  lying  traps  or  directly  from  the  valence 
band  (Zener  breakdown). 

(3)  The  electron  velocity  becomes  so  large  that 


electron-hole  pairs  can  be  generated  through  collisions.  These 
in  turn  are  again  accelerated  by  the  electric  field  so  that 
additional  electron-hole  pairs  are  created.  As  a  consequence 


s  single  carrier  csn  produce  sn  avalanche  of  electron- 
cole  pairs  {&7alscci;  tossaoxn). 

Junction  temMoxa  by  tssmsllsg  csn  be  expected  to  be 
one  of  the  chief  breahdGwn  mechaslsrs  In  the  reverse  biased 
xetal-sexiiccaaductor  diode.  Ibis  fact  be  cones  apparent  by 
sxBmimlng  Fig.  5.2  for  large  electric  fields  at  the  contact. 

I  Is  greatly  increased  dme  to  fff  and  breakdown  of  the 
j unction  may  occur . 

Junction  breakdown  by  the  Zener  effect  Is  less  likely 
to  occur  In  metal— semiconductor  devices  than  In  F-2£  junctions. 
Ibis  Is  due  to  the  fact  that  current  flow  In  the  metaX- 
semlccnductor  device  Is  chiefly.  If  not  exclusively,  by 
majority  carriers.  Since  /m~Xe  is  typically  much  less  than 
S_  (the  original  assumption  on  which  the  approximation  that 
Ck»o  Is  based)  tunneling  by  the  Zener  mechanism  Is  greatly 
reduced  for  a  majority  carrier  device. 

Junction  breakdown  by  the  avalanche  mechanism  Is  also 

of  prime  importance  for  a  metal- semiconductor  contact  under 

large  reverse  bias.  At  very  large  reverse  biases  it  may 

happen  that  an  electron  emitted  over  the  barrier  from  the 

metal  into  the  semiconductor  may  gain  more  energy  due  to 

the  high  electric  field  during  collisions  that  it  loses 

because  cf  a  collision.  Thus,  the  kinetic  energy  of  the 

electron  will  continue  to  increase  until  it  exceeds  E  • 

S 

At  this  point  a  collision  can  generate  electron-hole  pairs, 
which  in  turn  may  produce  additional  electron-hole  pairs 


and  the  avalanche  process  has  begun,  sender  this  condition 


the  reverse  current  I _  (reverse  current  ssrlor  to  cr&asdcwn) 

so 

is  greatly  increased  fcj  the  addition  of  free  carriers,  i.e.. 


T~  -  357  (5e5*i) 

o  so 

In  which  E  is  the  avalanche  nmitipii cation  factor.  At  a 
certain  fieM  strength,  which  is  constant  for  a  given  semi¬ 
conducting  material,  M  be  cones  infinite  and  breaMcxn  occurs, 
empirically  it  is  found  that  X  may  he  expressed  as 


(6.5.2) 


in  which  7,  is  the  reverse  voltaxe  at  which  bread-down 
or 

occurs  and  33  is  a  relatively  large  (5  to  6)  constant  dependent 
on  the  particular  seniconducting  material. 


6.6  Qualitative  Comparisons  Between  the  Theoretical 
I— V  Characteristic  and  Experimental  Measurements. 

The  complete  theoretical  expression  for  the  I-V  character 
istic  was  summarized  in  Section  6,4  for  an  idealized  metal- 
semiconductor  contact.  Although  the  complete  expression 
would  be  quite  difficult  to  evaluate,  it  could  conceivably 
be  done  and  result  in  a  theoretical  plot  of  the  I-V 
characteristic  by  assuming  the  necessary  constants  and 
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techni ones.  However,  without  aaSisg  a  direct  co-par i sen 
betxssn  a  theoretical  plot  of  this  type  ard  cal  experimentally 
determined.  characteristic,  it  is  possible  to  rase  certain 
qualitative  comparisons  between  experimental  measurements 
arid,  the  general,  fern  of  3q.  6.4.1.  5his  type  of  comparison 
will  be  undertaken  in  this  section  in  hopes  of  Justifying  at 
least  the  form  of  the  theoretical  expression. 

In  FIs:,  6.4  an  experimentally  determined  I-V  character¬ 
istic  Is  shown  for  a  Pt~5i  (2*-Type)  contact.  Bor  this 
particular  application  it  is  convenient  to  plot  2n{Ij.) 

iu 

versus  ¥  .  Using  Bo.  6.4,1,  InJI^)  becomes 


in(I,J  -  :n'l*exr(-eV<,/£T>]  -  ln[e^p( eV^ArD-1  ;{6.6.1 ) 


which  under  forward  bias  conditions  nay  be  written  as 


Ln(I. )  -  ln[l'ex?(-eV./KT)]  t  eV^/KT  (6.6.2) 

«  C  U  A 

since  exp(eVf/KT)»l.  Thus,  if  ln(It)  is  plotted  against 
a  straight  line  should  result  with  slope 


/•»»*> 
--  e/ri.i 


(6.6.3) 


Experimentally  it  is  found  that  for  this  particular  contact 


0 


1 


2 


3  .*■ 
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FIGURE  6.4: 


Current-Voltage  Characteristic  for  a  Ft-Si  (N-Type) 
Contact.36  Diara.=1.54(10"2)  era. 


In  which  5-1.02.  S>Is  saaa.ll  deviation  from  the  slops  as 
predicted  by  3a,  5,6*3  can  be  attributed  to  the  ncnideal 
nature  of  the  contact,  or  aaore  specifically  to  the  presence 
of  a  thin  lnter£aclal  layer  of  foreign  material  and/or 
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nommiforsa  contact  between  the  metal  and  semiconductor, 
in  general,  $  is  usually  included  in  the  expression  for  the 
I-V  characteristic  to  account  for  the  nonideal  nature  of  the 
contact  and  can  be  expected  to  be  slightly  different  for 
each  derice.  If  s  is  included,  it  say  be  determined  by 
the  following  relation: 


5  ~  KT  dln(l; 


( V_  sufficiently  large) 


(6.6.5) 


The  deviation  from  the  constant  slope  (m^.)  at  higher 
current  [ln(I.t)>10"’if  amps]  can  be  attributed  to  the  series 
resistance  associated  with  the  bulk  semiconductor#  Also, 
by  examining  Eq.  6,6.2  a  means  of  determining  I^exp(-eVd/KT) 
becomes  apparent*  This  parameter  may  be  determined  by  an 
extrapolation  of  the  ln(It)  versus  curve  until  inter¬ 
section  occurs  at  the  lndj.)  axis  (V^O)  as  shown  in 


The  exponential  dependence  of  I*exp(-eVd/KT)  on  the 


-  *v  -t  f 


-J, 


equilibrium  diffusion  potential  '7^)  can  be  shown  by  a  plot 

of  tbs  type  shown  in  Pig,  6*5.  Here  Inline  “  d/Ai)  is 

plotted  against  Vd  for  severs!  netsl-Si  (H-Type )  contacts, 

showing  a  constant  slope  of  appr oxi ra  te  ly  e/KT.  Although 

exact  agreement  Is  not  shown,  the  deviation  can  be  attributed 
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to  inaccuracies  in  the  measurement  of 

The  temperature  dependence  of  I^exp  ( -eV^/XT )  can  be 
shown  by  a  plot  of  the  type  shown  in  Fig,  6,6.  Here 
In ( I  )  versus  lfl  Is  plotted  under  forward  bias 

for  a  Au-Si  (M-Type)  contact  with  eVd=0.08.  The  slope  is 
constant  and  equal  to 

5a  «  -  e(Vd  -  7f)  =  0.799  ev  (6.6.6) 

E_  is  usually  referred  to  as  the  thermal  activation  energy 
for  the  contact. 

The  voltage  dependence  of  the  I-V  characteristic  at 
high  reverse  bias  is  shown  by  Fig.  6.7.  Here  ln(Xr)  versus 
V  is  plotted,  where 

Ir  =  Iorexp(yb)  (6.6.7) 

and  y£  is  given  by  Eq.  4.47.  Ip  is  interpreted  as  the 

effective  reverse  current  flowing  in  this  region  {large 

reverse  bias)  of  operation  and  J,or  is  the  effective  satura- 

42 

tlon  value  of  Ir#  For  silicon  €^=11.7€0  ,  and  good  agreement 


is  shown  between  the  theoretical  value  predicted  by  image 
effects  and  the  data  shown  in  Fig.  6.7. 


10"' 


Vd  (volts) 


.50  .55  .60  .65  .70  .75  .80  .85 

FIGURE  6,5*  Dependence  of  the  Saturation  Current  Density  on 

the  Equilibrium  Diffusion  Potential  (V,)  for 

a 

Several  Metal-3i  (N-Type)  Contacts. ^9 


i: 

L 


103/T  (°K) 


FIGURE  6.6:  Temperature  Dependence  of  the  Saturation  Current 
Un^er  Forward  Bias  for  a  Au-Si  (N-Type)  Contact f 
(eV.=0.08  ev  and  E  =0.79?  ev) 


CHAPTER  VII 


SUEKARY.  CONCLUSIONS  AND  RECOMMENDATIONS 
FOR  FURTHER  STUDY 

This  study  originated  as  an  effort  to  answer  the 
question  of  why  the  experimental  method  of  measuring  capaci¬ 
tance  as  a  function  of  reverse  bias  could  be  used  to 
accurately  predict  the  doping  concentration  and  yet  could 
not  be  expected  to  give  a  reliable  measure  of  the  eq'  •  librium 
diffusion  potential,  even  though  both  quantities  resulted 
from  an  interpretation  of  the  same  data.  Since  the  capaci¬ 
tance  technique  Is  based  on  the  Schottky  model  for  the 
contact,  the  most  logical  approach  to  explaining  this 
apparent  Inconsistency  was  to  examine  more  closely  the 
theoretical  expression  for  junction  capactiance,  or  more 
basically  the  Schottky  model  itself,  while  paying  particular 
attention  to  any  approximations  which  could  limit  accuracy. 
The  Schottky  model  is  based  on  the  assumption  that  the 
space  charge  associated  with  the  depletion  region  formed 
at  a  metal-semiconductor  contact  is  constant  and  due 
entirely  to  the  uniform  and  complete  ionization  of  donor 
Impurities  within  the  depletion  region.  Outside  the 
depletion  region  charge  compensation  due  to  the  presence  of 
free  electrons  is  complete  so  that  the  net  space  charge  is 
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zero*  A  core  accurate  model  for  space  charge  In  the 
depletion  region  is  developed  In  Chapter  IV  which  considers 
the  possibility  of  incomplete  and  nonuniform  donor  ioniza¬ 
tion  and  the  partial  compensation  to  positive  space  charge 
due  tc  the  presence  of  free  electrons  in  the.  immediate 
vicinity  of  the  depletion  region  edge.  Once  the  space 
charge  density  is  known  as  a  function  of  distance  in  the 
depletion  region,  Poisson’s  equation  (this  study  assumes  a 
one-dimensional  application)  can  be  used  to  solve  for 
expressions  for  the  electric  field  and  potential  variation 
within  this  region.  A  comparison  between  the  solutions 
based  on  both  models  shows  that  the  Schottky  model  is  quite 
accurate  in  the  region  near  the  contact;  however,  there  may 
be  significant  error  in  the  region  near  the  depletion  region 
edge.  Furthermore,  this  error  increases  as  the  degree  of 
ionization  decreases. 

Sven  though  there  may  be  large  differences  between  the 
electric  field  and  potential  plots  for  the  two  models  on  a 
point  by  point  basis,  a  comparison  between  the  two  space 
charge  density  plots  reveals  quite  clearly  why  the  Schottky 
model  can  give  accurate  results  for  certain  applications.  A 
comparison  of  this  type  shows  that  the  Schottky  model 
represents  a  fairly  accurate  estimate  for  the  effective 
space  charge  density  since  the  error  predicted  by  the  more 
exact  model  on  a  point  by  point  basis  is  compensating  and  the 
overall  difference  between  the  two  models  can  be  quite  small. 


if  the  degree  of  ionization  is  fairly  close  to  unity.  This 
is  equivalent  to  saying  that  the  areas  under  the  two  space 
charge  density  plots  (Flg^s  4.2b  and  4.3b)  are  approximately 
the  same. 

In  Chapter  V  the  differential  capacitance  technique  for 
measuring  the  equilibrium  diffusion  potential  was  examined 
in  detail*  Using  the  expressions  for  Junction  capacitance 
(C^“^)  based  on  both  models  the  diffusion  potential  and 
Intercept  voltage  are  related  as  follows: 

7d  =  Vo  (Schottky  Model)  (7.1) 


V,  =  V  +  —(y-lny)  (More  Exact  Model)  (7*2) 

Comparing  these  two  expressions  one  can  expect  the  diffusion 
potential  predicted  by  Eq.  7*1  to  be  in  error  by  at  least 
KT/e  (using  Eq.  7.2  with  y=l )  and  this  error  will  increase 
as  the  degree  of  ionization  decreases.  For  y=0.05  the 
error  is  approximately  3(KT/e);  however,  this  is  not  a  true 
representation  of  the  actual  error  since  the  temperature  must 
also  decrease  as  y  decreases. 

Another  interpretation  of  the  data  from  a  capacitance 

_2 

(Cj  )  versus  reverse  bias  voltage  plot  results  in  an 
expression  for  the  doping  concentration.  In  Chapter  V 
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Nd  was  shown  to  be  proportional  to  the  slope  and  the 
predicted  slopes  from  both  models  were  Identical  fcr  suf¬ 
ficient  reverse  bias  (yQ>4).  The  theoretical  expression  for 
Nd  for  both  models  Is  given  by 


N.  - - 

a  e£  d(C“2) 


(yc>4) 


(7.3) 


and  is  known  to  agree  quite  well  with  results  obtained  by 


other  means  of  measuring  N^.  Why  the  apparent  inconsistency 
between  i?d  and  Vd  using  the  differential  capacitance  method 
based  on  the  Schottky  model?  The  inaccurate  measurement  of 


Vd  by  interpreting  Vd  as  equal  to  VQ  is  not  the  fault  of  the 
model,  but  results  from  a  false  interpretation  of  VQ.  The 
experimental  method  of  determining  Vq  is  by  extrapolating  a 
Cj"  versus  Vr  plot,  which  is  linear  as  shown  in  Eq.  7«3» 
to  a  point  where  it  intersects  the  Vr  axis,  i.e.,  VQ.  This 
straight  line  passes  through  a  region  of  small  reverse 


bias  and  the  eventual  value  of  VQ  falls  in  the  region  of 
negative  Vp  (see  Fig.  5.4a).  The  apparent  assumption  is  that 
the  capacitance  (C^*"2)  remains  linear  in  this  region  and 
this  is  contrary  to  fact  as  shown  in  Fig.  5.5.  The  non- 
linear  variation  in  capacitance  (C^  )  at  small  values  of 
reverse  bias  can  be  explained  quite  easily  by  examining  the 


term  corresponding  to  the  free  electron  concentration  in 
Eq.  4.23,  or 


( 0<x<W ) 


(?.*) 


n(x)  =  n  exp(-y ) 

G 


Under  a  fairly  large  reverse  bias  yQ  is  large  and  y  is 
increased  so  that  n(x)«0,  which  is  consistent  with  the 
assumption  on  which  the  Schottky  model  is  based;  however, 
as  V_  is  decreased,  yQ  becomes  smaller  and  eventually  a 
point  is  reached  at  which  one  may  no  longer  neglect  the 
free  electron  population  at  the  depletion  region  edge. 

This  value  of  yQ  can  be  Interpreted  from  the  capacitance 

-2 

(Cj  )  versus  voltage  plot  as  shown  in  Fig,  5,5,  Since  the 
capacitance  (C^"2)  becomes  nonlinear  for  all  values  of  yQ 
less  than  four  (V^-s-v^-0,104  volts)  and  all  values  of  y, 
then  one  may  interpret  this  as  the  point  at  which  the 
Schottky  model  assumption  that  the  depletion  region  is  free 
of  electrons  is  no  longer  valid.  There  is  a  small  variation 
in  this  point  due  to  y,  which  is  consistent  with  this  line 
of  reasoning  since  n„  in  Eq,  is  a  function  of  y,  i,e,, 
nc=YNd*  Thus,  if  y  decreases  the  free  electron  population 
at  the  depletion  region  edge  will  also  decrease  and  the 
cut  off  value  of  yQ  can  be  expected  to  increase.  This  is 
shown  quite  clearly  in  Fig,  5,5, 

Since  the  interpretation  of  VQ  is  based  on  the  intercept 
value  of  a  linear  plot,  the  value  of  Vo  may  be  in  error  and 
the  plot  can  still  predict  the  correct  slope.  This  is 
pointed  out  by  Fig.  5,5  in  which  the  Schottky  and  more  exact 

a 
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models  predict  the  same  slope  for  yQ>4,  but  there  Is  signif¬ 
icant  difference  in  the  values  predicted  for  V  .  Furthermore, 
the  slope  is  Independent  of  y  (yo>ii)  and  the  intercept  value 

of  V  is  not.  The  reason  is  again  due  to  the  correct  Inter- 
o 

pretation  of  the  components  of  space  charge  for  different 
degrees  of  lon'zation.  If  all  donor  atoms  are  assumed  to  be 
ionized  and  yo>^»  then  the  slopes  of  the  Schottky  and  more 
exact  model  are  in  complete  agreement;  however,  if  yd  then 
the  positive  space  charge  furnished  by  positive  donor  sites 
is  reduced  and  cne  would  expect  an  increasing  disagreement 
between  the  two  models  for  space  charge  density  as  y  decreases 
This  is  shown  in  Fig.'s  k.2b  and  4.3b.  Furthermore,  since 
the  effective  depletion  width  is  reduced  for  smaller  values 
of  v  one  would  expect  larger  values  of  capacitance  (assuming 
a  constant  value  of  bias)  for  lower  values  of  y  (see  Fig,  5,5) 
Since  a  lower  value  of  y  implies  a  larger  value  of  capacitance 
for  the  same  bias  voltage,  one  would  also  expect  the  intercept 
value  of  V  (based  on  a  linear  estrapolation)  to  be  y  depen- 
dc’'t,  with  V  increasing  as  y  decreases.  This  is  also  evident 

O  w 

rroo  Fig.  5,5. 

The  agreement  between  the  two  models  in  predicting  Nd 
for  yQ>^  stems  from  the  fact  that  the  two  models  are  com¬ 
pletely  equivalent  in  this  range  of  bias  with  regard  to 
capacitance.  The  voltage  dependence  of  the  junction  capaci¬ 
tance  arises  from  the  change  in  stored  charge  within  the 
depletion  region  with  and  incremental  change  (C^=dQ/dV)  in  the 


applied  bias  and  this  addition  and  subtraction  of  charge 
will  be  independent  of  the  initial  stored  charge  (pro¬ 
portional  to  y)  as  long  as  the  depletion  region  is  free  of 
charge  compensating  electrons.  For  the  regior  of  bias 
voltage  in  which  the  free  electron  population  at  the 
depletion  region  edge  can  no  longer  be  neglected  (yQ<4),  the 
change  in  capacitance  associated  with  a  change  in 

bias  is  no  longer  constant  since  the  density  of  compensat¬ 
ing  electrons  varies  exponentially  (Eq.  7 • 3 ) •  Furthermore, 
the  electron  population  is  dependent  on  y  so  that  the  slope 
can  be  expected  to  be  altered  by  an  additive  constant  due  to 

the  influence  of  y.  In  addition,  the  change  in  capacitance 

—2 

{ C j  )  caused  by  a  change  in  y  Will  eventually  reverse 

Itself  in  this  region,  i.e.,  higher  values  of  y  will 

produce  lower  values  of  capacitance  (assuming  the  same  value 

of  bias)  since  the  electron  population  at  low  values  of  y 

is  more  strongly  influenced  by  y  than  by  y  (y  approaches  a 

minimum  value  which  is  proportional  to  V^).  The  above 

conclusions  are  illustrated  quite  well  in  Fig.  5.5. 

In  conclusion,  the  inaccuracy  introduced  by  using 

Eq.  7.1  to  interpret  the  diffusion  potential  from  a 
••2 

capacitance  (C j*"  )  versus  voltage  plot  is  due  to  the  mis¬ 
representation  of  the  Schottky  model.  The  Schottky  model 
assumes  that  the  depletion  region  is  free  of  electrons  and 


yQ  has  a  minimum  value  for  compliance  with  this  assumption. 

—2 

Since  capacitance  (Cj  )  is  a  nonlinear  function  of  bias 
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below  this  minimum,  a  linear  extrapolation  for  VQ  produces 
significant  error  in  V^.  However,  if  Eq.  ?.2  is  used  to 
interpret  VQ,  then  one  may  expect  an  accurate  value  for  the 
equilibrium  diffusion  potential. 

The  influence  of  image  effects  on  the  accuracy  of 
Eq.  7.2  was  also  considered  in  Chapter  V  and  the  conclusion 
drawn  that  image  effects  may  be  safely  neglected  for  this 
particular  application.  This  is  due  to  the  fact  that  the 
reduction  in  barrier  height  caused  by  image  effects  is  quite 
small  in  comparison  to  the  barrier  height  (neglecting  image 
effects)  at  the  contact,  i.e.,  yQ~y {f*y0  a t  all  values  of 
bias  even  though  y£  increases  with  bias. 

This  study  has  also  considered  a  theoretical  develop¬ 
ment  of  the  current-voltage  characteristic  of  an  idealized 
contact  based  on  a  diffusion  model.  Although  justification 
for  the  accuracy  of  the  I-V  model  developed  on  a  quantitative 
basis  has  not  been  attempted,  a  qualitative  comparison 
between  the  model  and  experimental  measurements  have 
produced  favorable  results.  The  model  is  able  to  predict 
the  correct  slope  for  the  forward  characteristic  within  a 
small  multiplicative  constant  and  this  deviation  can  be 
attributed  to  the  nonideal  nature  of  the  contact.  Further¬ 
more,  the  forward  characteristic  shows  the  correct  temperature 
dependence  and  the  reverse  saturation  current  agrees  well 
with  the  experimentally  observed  exponential  dependence  on 
the  equilibrium  diffusion  potential.  Also,  the  model  is 
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able  to  predict  a  voltage  dependence  of  the  reverse 
characteristic  which  agrees  well  with  experiment.  This 
dependence  arises  from  the  increasing  importance  of  image 
effects  at  large  reverse  biases.  The  theoretical  model 
predicts  a  high  probability  of  breakdown  due  to  tunneling 
which  is  also  consistent  with  experimental  results.  In 
short,  although  the  theoretical  I-V  model  would  require 
fairly  sophisticated  numerical  techniques  for  evaluation, 
it  should  agree  quite  well  with  experimental  results  on  a 
quantitative  basis. 

Since  the  results  of  this  study  are  quite  conclusive 
with  respect  to  the  use  of  the  differential  capacitance 
technique  for  determining  the  diffusion  potential  and 
doping  concentration,  this  method  should  compare  favorably 
with  values  of  measured  by  the  photoresponse  technique. 
Thus,  an  experimental  study  of  this  type  would  warrant 
consideration.  Also,  it  would  be  interesting,  although 
the  practicality  may  be  questioned,  to  pursue  a  complete 
numerical  evaluation  of  the  theoretical  current-voltage 
characteristic  for  the  purpose  of  comparing  the  results  with 
experimental  measurements. 
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APPENDIX  A 

THE  KOTT  MODEL  OP  A  KETAL-SSKIC  GNDOCTOH 
RECTIFYING  CONTACT 

The  Mott  Eodel^  for  a  metal-  (:.r-?ype  ^semiconductor 
contact  will  not  be  discussed  In  as  such  detail  as  the 
Schottky  and  Bethe  node  Is  since  Its  application  is  limited 
to  only  a  very  thin  barrier.-  Furthermore ,  it  represents  a 
special  case  of  the  Schottky  model  which  Is  discussed  in 
Chapter  III.  The  model  and  equations  based  on  the  model 
are  Included  only  for  completeness,  although  Its  Importance 
should  not  be  overlooked  since  It  was  one  of  the  first 
successful  models  which  could  predict  and  mathematically 
account  for  the  rectification  phenomena  observed  at  a  metal- 
semiconductor  contact. 

The  potential  distribution  for  the  Mott  barrier  is 
shown  in  Fig.  A,t.  The  barrier  Is  defined  as  a  terrier 
which  extends  throughout  the  U-Type  semiconductor,  or  In 
other  words  .joins  both  metal  contacts  of  the  device.  The 
left-hand  contact  is  regarded  as  ohmic  since  the  terrier 
at  x=W  offers  very  little  opposition  to  current  flow.  The 
basic  assumption  of  the  Kott  model  is  that  due  to  the  nature 
of  the  barrier  too  few  Impurities  of  the  semiconductor  are 
ionized  to  disturb  the  electric  field  in  the  semiconductor. 


3 


as 

5 
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Since  the  charge  associated  with  the  system  must  be  assumed 
by  the  metal  surfaces,  the  electric  field  is  constant  in  the 
semiconductor  and  the  potential  function  is  thus  linear. 
Stated  mathematically,  Poisson *s  equation  takes  the  fora 


(G^x=W) 


(A.l) 


since  there  is  no  space  charge  associated  with  the  semi-¬ 
conductor.  Upon  integration  and  evaluation  at  the  boundaries 
one  obtains  an  expression  for  the  electric  field,  or 


dV  _ 
dr 


-  six) 


v 


(A.2) 


A  second  integration  yields  an  expression  for  the  potential 
distribution  of  the  semiconductor,  or 


V(x)  = 


V ,  -  V 
d  _ a 

’rf 


+  C 


(A.3) 


in  which  C  say  be  evaluated  by  noting  that 


i 


vio )  =  -E-D/e  =  -Vb  (A.4> 

Thus,  the  potential  is  linear  and  may  be  expressed  as 


C 
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TEE  3ETHE  KODEL  OF  A  KSTAL- SEMICONDUCTOR 
RECTIFYING  CONTACT 

The  Bethe  model^  for  a  netal-semiccnductor 
contact  is  shown  in  Fits*  3.1.  The  assumptions  used  to  model 
the  contact  are  identical  to  those  of  the  Schottky  model 
with  the  additional  assumption  that  the  interface  is  an 
insulating  layer  of  thickness  r,  separating  the  metal  from 
the  semiconductor#  The  model  is  particularly  applicable  to 
a  semiconducting  Gxide  layer,  which  may  have  been  deliberately 
applied  to  the  semiconductor  before  forming  the  Eetal  contact, 
or  may  have  resulted  from  the  process  of  applying  the  metal 
contact.  In  either  case  the  semiconducting  oxide  forms  a 
thin  insulating  layer  and  the  model  assumes  the  transition 
from  the  oxide  to  the  pure  semiconductor  Is  abrupt# 

Under  the  above  assumptions  Poisson1 s  equation  takes 
the  following  form  for  a  planar  contact: 


daV(x)  _  0 
dx2  ' 


(0<x<r) 


(B.l ) 


for  the  region  of  the  insulator,  and 


dxa 


(r<x<W) 


(B„2) 


for  the  region  of  the  pure  semiconductor  with  completely- 
ionized  donor  concentration  The  potential  within  the 

barrier  may  be  found  by  solution  of  Poisson* s  equation, 
subject  to  the  following  boundary  conditions: 


fi)  V(W)  =  vd  +  vr 


CB.3) 


(2)  V(0)  =  o 


(B.4) 


(3)  =0  at  x=W 


(B.5) 


with  the  potential  and  electric  field  continuous  at  x=r. 
The  solution  may  be  written  as 


V(x)  =  — -r)x  (0<x<r) 
s 


(B.6) 


V(x)  -  pp-[W2-rs-(W-x)3]  (r<x<W) 


(B#7 ) 


The  Be the  model  gives  the  following  result  for  the  barrier 


4  -  ,  *  C*_ 


width  W: 


I  r 

,  i 

I  r 

o, 

■i  I 


w*  - +  a“<W 


(B.8) 


Since  the  Bethe  model  neglects  electrons  and  holes  as  part 
of  the  space  charge  of  the  barrier,  the  charge  per  unit  area 
is  given  by 


Q(W+r)  =  eNd(W~r) 


(B.9) 


giving  the  following  result  for  barrier  capacitance s 


"  W  - 


2C  -I? 

r2  +  a#  W 


(B.10) 


One  should  note  that  the  conditions  of  the  problem  as  stated 
in  Poisson* s  equation  make  the  solutions  physically  valid 
only  for  values  of  the  barrier  width  equal  to,  or  greater 
than  the  thickness  of  the  insulating  layer. 

As  a  means  of  comparing  the  Bethe  and  Schottky  models 
the  junction  capacitance  versus  reverse  bias  voltage  is 


shown  in 


£,2«  The  extrapolation  indicated  in 


Illustrates  how  one  may  theoretically  estimate  the  thickness 
of  the  insulating  layer  through  the  intercept  on  the  voltage 
axis.  For  an  applied  forward  bias  voltage  exceeding  the 
equilibrium  diffusion  potential  (V^),  the  barrier  should 
behave  as  an  ordinary  capacitance  with  a  dielectric  layer 
of  thickness  r.  This  would  correspond  to  the  region  of  the 


a-LihiLi 


FIGURE  fl. 2t  Comparison  of  the  Bethe  and  Schottky  Model 

-2 

Plots  of  C  versus  V^« 


horizontal  dashed  line  in 


However,  the  diffusion 


Fig,  B.2. 

of  charge  carriers  into  and  beyond  the  barrier  does  not 
attain  a  steady  state  condition  until  a  certain  time  has 
elapsed  and  this  time  may  be  on  the  order  of  time  needed  to 
charge  the  barrier.  This  would  have  the  effect  of  either 
causing  the  measured  capacitance  to  fall  below  or  above  the 
constant  value  Indicated  in  Fig.  B.2. 
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APPENDIX  C 


AN  APPROXIMATION  FOR  THE  QMTC  OF  A 
PARABOLIC  POTENTIAL  BARRIER 

Kemble* ^  has  approximated  the  transmission  coefficient 
(QMTC)  for  electrons  of  momentum  s  incident  on  a  parabolic 
potential  barrier  with  the  presence  of  an  image  force  and 
with  reference  to  Fig.  C.l  his  results  are  as  follows: 


Case  I:  QMTC  = 


(C.l) 


Case  II:  QMTC  ,  —,4^) 


(E'>g  ) 
n 


(C  *2 ) 


Case  III:  QMTC  =  | 


(E'=3  ) 
r.  m 


(c.3) 


where 


s  *En 


:(;<)  =  ^2m*[E;rP(x)]J2 


(C.4) 


(C.5) 


r,.  •  *•_  v 
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ex?(E,/ET) 


(HcAd)(1-Y) 


2y= 


0.1) 


Equation  D.l  can  be  arranged  to  fora  a  quadratic  equation 
In  y,  or 


ay2  by  -  c  =  0 


0.2) 


where 


a  -  2exp(E-/K?) 


0.3) 


and 


b  = 


=  Jf  /*2 

c/:d 


(D.4) 


Using  the  standard  quadratic  formula,  the  principal  value 
of  y  say  be  expressed  as  follows: 

(-1  +  [1  +  8{N  /K,)exp{2./KT)]) 

y  =  £ - = - U-  (D.5) 

^N^expCS^/KT) 


The  effective  density  of  states  for  the  conduction  band 
(Nc)  may  be  expressed  as^ 


N 


c 


_  3/2 

4/2(mn*KT) 

e 

- 3 - 

h 


(D.6) 


5GK: 
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( F,  As,  and  Sb)  to  produce  K-Type  silicon  and  gerraniun, 

essentially  all  inpurities  will  be  ionized.  This  is  due  to 

the  relatively  low  ionization  energies  associated  with  these 

44 

dopants,  ranging  from  0.001  to  0.055  e7»  However,  if  one 

I>c 

considers  B-Type  SiC  using  A1  doping,  S, js0.2;>  sv  ■*'  and 
ionization  of  impurities  is  relatively  incomplete  at  room 
temperature. 

(2)  The  doping  density  (Nd)  and.  the  effective  electron 
mass  (n*)  also  influence  the  character  of  y»  although  this 
influence  nay  gc  unnoticed  because  of  the  importance  of 
and  T.  However,  one  nay  draw  the  general  conclusion  that 
Y  Is  proportional  to  the  effective  mass  and  inversely 
proportional  to  the  doping  concentration. 
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APPENDIX  E 

VERIFICATION  OF  EQUATION  4.32  AS  TEE 
SOLUTION  OF  EQUATION  4.27 

In  order  to  show  that  Sq.  4.32  is  a  solution  of  Eqc 
4.27  It  must  be  shown  that  Ea.  4.32  satisfies  the  differ¬ 
ential  equation  (Sq.  4.27)  and  meets  the  boundary  condition 
of  Sq.  4.30.  This  result  can  be  shown  quite  easily  by 
performing  the  following  differentiations: 


Upon  substitution  of  E.2  and  E.3  into  E.l  one  obtains 


i  i 

r 


a2y  _ 


1  +  (l/Y-l)e"y 


-  ye-y 


(B.4) 


which  is  identical  to  Eq.  4# 2?.  Thus,  Eq.  4,32  satisfies 
the  differential  equation, 

Next,  considering  the  boundary  condition: 


*£  - 


=  0  at  y=0 


f^|y=0  =  -^2[ln(Ye°+l  -  Y)  +  ye°  -  y]^ 


(2.5) 


(2.6) 


|y=0  =  -/2[ln(l )  +  y  -  y]^  = 


(2.7) 


so  that  Eq.  4.32  satisfies  the  boundary  condition  of  Eq. 
4.30.  Thus,  since  Eq.  4.32  meets  both  requirements  as 
stated  above  it  must  be  considered  a  valid  solution  to 
Eq.  4.27  subject  to  the  boundary  condition  of  Eq.  4.30. 
q  »e ,d» 
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APPENDIX  F 

THE  RUNGE-KUTTA  METHOD  FOR 
NUMERICAL  INTEGRATION 

The  solution  of  Eq.  4.32  is  accomplished  by  a  numerical 
integration  employing  the  Runge-Kutta  technique.  In  general, 
this  technique  takes  the  interval  (en, 8n+i }  and  breaks  it 
into  two  or  more  subintervals.  The  integral  of  the  function 
f (a , a)  over  the  whole  interval  is  then  calculated  as  the 
sum  of  the  Integrals  over  the  sub intervals.  The  function  is 
taken  to  be  constant  over  each  subinterval  and  by  Judicious 
choice  of  the  points  at  which  the  function  is  evaluated  a 
low  truncation  error  can  be  obtained.  In  addition,  when 
using  the  Runge-Kutta  method  only  one  Initial  point  (®0»A0) 
is  needed  for  the  iteration  process  to  begin.  The  chief 
disadvantage  of  this  method  is  that  it  requires  several 
evaluations  of  f($,A)  for  each  point  of  integration  which 
makes  it  somewhat  slower  than  other  methods. 

The  specific  Runge-Kutta  method  used  for  the  solution 
of  this  problem  is  of  order  four,  or  the  interval  (n)  is 
divided  into  four  subintervals.  The  integration  is  then 
given  by  the  following  system  of  equations: 


Z0  =  nf(«n,4n) 


(P.l) 


r.  =  nf(e  + 


n 


aIfl,An 


P10Z0} 


(F»2) 


=  n  r(e  -> 

z  ^ 

*•  •< 


a2n’  An 


P  T 
H20*0 


l21Zl 5 


(F.3) 


i0  -  nr(8 


n 


'  a3n,in 


‘Vo  +  e31rl  +  S32I2>  (P-4) 


ArHl  "  An  '*■  (aZ0  +  fcZl  *  cZ?  +  dI3^  (F.5) 

Although  a  large  number  of  parameters  appear  in  Eq.’s  F.l 

through  F.5  they  may  be  determined,  at  least  in  part,  by 

equating  Eq,  F.5  to  the  Taylor-serles  expansion  of  ^n+1 

about  the  point  ( 9  ,a  )•  This  can  be  done  so  as  to  achieve 

n  n 

4 

agreement  through  terms  in  n  ,  yielding  solutions  that  have 
an  error  of  approximately  (n^). 


9 


1 


One  possible  solution  is 
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z0  =  nf(  W 


(P.6) 


Z1  11  fif(®n  +  Q/2,An  +  n/2) 


(P.7) 


s  nf  (•  +  n/2,A  +  n/2) 

w  n  n 


(P.8) 


x~  =  nf(<?  n 


n 


n 


x2) 


(P.9) 


R+i  =  %  +  1/6(I0  +  2I1  -  2IS  +  V 


(P.10) 


where  dA/d«=f (e , a)  with  the  initial  starting  point  specified 
as  (®0»A0).  Here  the  interval  n  is  divided  into  sub¬ 
intervals  of  n/6,  n/3»  n/3*  and  Cl/6,  The  function  is 
evaluated  at  the  left-hand  side  first,  then  twice  at 
extrapolated  center  points,  and  finally  at  the  extrapolated 
right-hand  side. 

The  complete  computer  program  used  to  integrate  Eq.  4,32 
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by  the  Runge-Kutta  method  of  order  four  follows.  In 
addition  to  calculating  the  normalized  electron  potential, 
the  program  is  designed  to  calculate  the  space  charge 
density  based  on  both  models  and  present  the  results  in  a 
normalized  form  for  comparison  purposes.  Also  Incorporated 
into  the  program  is  a  plot  subroutine  which  allows  the  data 
generated  to  appear  In  a  more  useful  form.  The  program  is 
written  in  Fortran-IV  for  use  in  the  GE-225  computer  and 
takes  approximately  three  minutes  running  time  for  each 
parameter  change. 
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SUBROUTINE  PLOTT  (N,N?,X,Y1,22,Y3) 

DIMENSION  X(100)  fYl(iOO)  ,Y2(100)  ,Y3(100) 

DIMENSION  A (100)  ,?LOT{100) 

COXXON  G  ,  BUiK  ,  PHD ,  STAR  ,EQ,  PLUS  .COM,  EXE 

7  F0BHA?(lX,2210.3,4Xf 95^1 ) 

8  FCBJSATdX^SlOO.^X.eSAl) 

9  POHHAT  ( IX , 4E1 0 . 3 , 4X , 75A1 ) 

300  P05XAT ( 1H1 . 45X ,  1 6HPLOT  SUBROUTINE  ) 

310  FOH*AT(lH  ,22X, 13,433  VALUES  0?  Di/DU(*},  V{=}5  K(+)  A3E 

1  PLOTTED,  83  AGAINST, 13, 313  VALUES  O?  U( . )  SCALE  FACTOR 

2  =,S10«3) 

KB=N? 

5>o  5  I=i,95 
5  ?LOT(  I  )=3L?JX 
XK=0.Q 
XXU=0.0 
DO  10  1=1, NP 
10  A(I)=Y1(I) 

KSW=2 
GO  TO  60 
40  DO  50  1=1, NP 
50  A(I)=Y2(I) 

M3W=3 

60  DO  90  1=1, NP 

IF(A(I)-XM)  70 ,75*75 
70  IF(A(X)-XMI)  80,90,90 
80  XMI=A(I) 


11  1 


GO  TO  90 
75  XK=ACI) 

90  CONTINUE 

IP(MSW-{N-1) }  20,^0,100 
100  NEA-94-(N-l )*1G 
3A=NRA 

SP=(XK-XMI)/RA 

R£P=A3S?{XKI/(XM-XMI ) )*RA+1 . 0 

NR=BEP 

PRINT  300 

PRINT  310,NP,NF,SF 
GO  TO  (320,330,340) ,N 
320  PRINT  325 

325  FORMAT* 7Xf4HX(.),5X,5BY!(*)) 

GO  TO  4 

330  PRINT  335 

335  FORMAT(7Xf4HX(.)f5X,5HYl(*),5X,5HY2(=)) 

GO  TO  4 
340  PRINT  345 

345  F0RMAT(7X,4HX(.),5X,5HY1(*),5X,5HY2(=),5X,5HY3(+)) 
4  PRINT  350 


350  FORMAT* IX, 116H, 


DO  290  1=1, NP 

PLOT(NR)=PRD 

L=Y1(I)/SF+REF 


,*>3 
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PLOT ( L ) =STAB 
GO  TO  (140,110,110) fN 
110  X=Y2 ( I )/SP+HEF 
PLOT(K)=EQ 
IP(N-2)  150,150,120 
120  J=Y 3(1)/ SF+EEP 
PLOT(J)=PLUS 
GO  TO  160 

140  PRINT  7,X(I),Yl(I),(PLOT(J)tJ=l,95) 
PLOT(L)=BLNK 
GO  TO  290 

150  PRINT  8,X(I),Y1(I),Y2(I),(PL0T(J),J=1,85) 
PLQT(L)=ELNK 
PLOT(K)=BLNK 
GO  TO  290 

160  PRINT  9, X(I)fYl(I),Y2(I),Y3(I),(PLOT(J),J=l, 75) 
DO  165  NO=l,95 
I65  PLOT ( NO ) =BLNK 
290  CONTINUE 
RETURN 
END 
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SUBROUTINE  FNCT  (X,YfF) 

COMMON  G , BLNK , PRD, STAR f EQ, PLUS , COM , EXE 
AA=SQRTF{ 2.0} 

BB=LOGF(G*EXPF( Y )+l . 0-G )+G*EXPF( -Y ) -G 
CC=SQRTF(BB) 

F=-AA*CC 


RETURN 


-1 

\ 


4 
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C  **********  RUNGE-KUTTA  METHOD  FOR  SOLUTION  OF  D.Ec  USING  GE415 

C  **********  E.  ALLEN  LINDSEY,  NOVEMBER  1968 

C  **********  THE  FOLLOWING  VARIABLES  ARE  USEDt 
****  A  AND  B  ARE  THE  INTERVAL  LIMITS 

****  h  IS  THE  STEP  SIZE 

****  y  IS  THE  INITIAL  VALUE  OF  THE  SOLUTION 
****  DELTA  IS  THE  INCREMENT  AT  WHICH  THE  SOLUTION  IS 
IS  TO  BE  PRINTED  OUT 
****  g  IS  THE  DEGREE  OF  IONIZATION 
DIMENSION  NEQ(80),P(1) 

DIMENSION  Xl{200) ,Y1(200) ,Y2(200) ,Y3(200) 

DIMENSION  X11(70)9Y11(?0),Y22(70),Y33(70) 

COMMON  G , BLNK , PRD, STAR, EQ, PLUS , COM , EXE 
999  FORMAT  (80A1) 

998  FORMAT  (//5X.19E  ON  THE  INTERVAL  U=,F8.4,6H  TO  U=,F8.4, 

1  17H  WITH  STEP  SIZE  =,  F 9.7/) 

997  FORMAT  (5X.51H  THE  TOTAL  NUMBER  OF  INTERVALS  FOR  THIS 
PROBLEM  IS  , 

1  16/) 

99 6  FORMAT  (5X,48H  THE  DEGREE  OF  IONIZATION  FOR  THIS 
PROBLEM  IS  G=, 

1  F8 . 4/ ) 

995  FORMAT  (5X.48H  THE  EXPECTED  ERROR  FOR  THIS  PROBLEM 
IS  E=(+OR-) , 

1  E16.8/) 

994  FORMAT  (8X,2H  Uf13X,2H  Y,15X,2H  V//) 

993  FORMAT  (5X,F12.7,3X,F13.8) 


3 


3 

1 
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992 

991 

990 

989 

988 


FORMAT  {5X,F12.?,3X,F13.8*3X,F13.8) 

FORMAT  (6F10.3) 

FORMAT  (//////) 

FORMAT  (2X,lAl,2X,lAl,2X,lAl,2X,lAl,2XflAi) 

FORMAT  (5Xt37H  ^NORMALIZED  DISTANCE  FROM  JUNCTION 

1  5X.33H  Y=N0RMALIZED  ELECTRON  POTENTIAL  / 

2  5X.35H  V=N ORMALIZED  SPACE  CHARGE  DENSITY  / 

3  5X.32H  W=SCHOTTKY  APPROXIMATION  FOR  V  ////) 

READ  989. PRD. STAR, EQ, PLUS, RLNK 


/ 


Nl=3 

c  **********  MAIN  PROGRAM 


1  READ  991,A,BtDELTA,HtY,G 
C  **********  CHECK  FOR  EOF 

IF  (H)  2,100,2 

2  NH=(B-A}/H 
PRINT  990 

READ  999 . ( NEQ( I ) , 1=1 , 80 ) 
E=H**5.0 


P(1)=SQRTF(2.0*Y) 
**********  PRINT  HEADINGS 
PRINT  999, NEQ 
PRINT  998,A,B,H 
PRINT  996, G 
PRINT  995, E 
PRINT  988 
PRINT  99^ 

PRINT  993, A,Y 


^c'aW«£S!j 
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c  **********  INITIALIZATION 
10  X=A 
1=1 
J=0 

c  **********  COMPUTATION  OF  THE  NUMBER  OF  PRINT  INTERVALS 
N=DELTA/H 
NN=N+1 
15  YO=Y 

c  **********  EVALUATION  OF  THE  FUNCTION  FOR  RUNGE-KUTTA  FORMULA 
30  CALL  FNCT  (X,YO,F) 

ZK1=F 
XH=X+H/2.0 
YO=Y+H*ZK1/2.0 
CALL  FNCT  (XH,YO}F) 

ZK2=F 

YO=Y+H*ZK2/2 , 0 
CALL  FNCT  (XH,YO,F) 

ZK3=F 

XH=X+H 

YO=Y+H*ZK3 

CALL  FNCT  (XH,YO,F) 

ZK4=F 

50  Y=Y+H* ( ZK1+2 . 0* ( ZK2+ZK3 )+ZK4 ) /6 . 0 
D=1 . 0+( 1 . O/G-l . 0 )*EXPF( ~Y ) 

V=  140/D-G*EXPF(-Y) 

C  **********  INCREMENT  ARGUMENT  AND  COUNTERS 
Z=I 


X=A+Z*H 

1=1+1 

C  **********  TEST  FOR  PRINT  INCREMENT 
I?  ( I-NN )  30,70(70 
70  PRINT  992, X, Y, V 
J=J+1 
Y1(J)=Y 
Y2(J)=V*Y1(1) 

IF  (X-P(l))  76,76,78 
76  Y3(J)=Y1(1) 

GO  TO  79 

78  Y3(J)=0.0 

79  X1(J)=X 

c  **********  TEST  FOR  FINAL  COMPUTATION 
IF  (NN-(NH+1) )  80,90,90 

80  NN=N+NN 
GO  TO  15 

90  S=(l. 0/DELTA )*B 
MN=S 

IF  (S-58.0)  96,96,92 
92  SCALE=S/58. 0+1.0 
NM=SCALE 
L=MN/NM 
K=-2*NM 
DO  94  JJ=1 ,L, 1 
K=K+NM 

Yll ( JJ )=Y1 (K+NM+1 ) 
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Y22(JJ)=Y2(K+NM+l) 

X33(JJ)=Y3(K+NM+l) 

94  Xll ( JJ)=X1 (K+NM-f-1 ) 

NP=L 
GO  TO  98 

9 6  CONTINUE 

DO  9?  J1=1,MN,1 
Y11(J1)=Y1(J1) 

Y22(J1)=Y2(J1) 

Y33(J1)=X3(J1) 

97  X11(J1)=X1(J1) 

NP=MN 

98  CALL  PLOTT  (N1  ,NP,X11 , Yll .Y22.Y33) 

C-O  TO  1 

100  CALL  EXIT 


END 


APPENDIX  G 


A  TYPICAL  COMPUTER  SOLUTION 


***************** 

DEPARTMENT . HE  * 

USER  NAME . .  ALl.ENL  INOSEY  * 

USER  MUmHER....  * 

■**•»**★*•<►■********* 


-  * 


*FJZ 


OCT  9  /OCT  9  ,  CARD  *AAU*  FLOATING  POINT 
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SUBROUTINE  PLOTT(N,NP>X, Y1,Y2»V3) 

DIMENSION  X{l00],Yl[l00),Y2[i00],Y3tlQ0] 

DIMENSION  AllOO!,  PLOT  £  10  0  3 

COMMON  G,BLNK,PRD*StaR>EO> PLUS, COM. EXE 

7  FORMaT[1X,2E10.3,4X,95Ai} 

8  FORMAT(lX,3Ei0.3,4X,85Aij' 

9  FORMATtlX,4E10.3,4X,75AU 

300  FORMaT [1H1.45X.16H PLOT  SUBROUTINE  1 

310  FORMaTIIH  ,22X,I3,43M  VALUES  OF  DY/DUt*),  Vt«J,  Wl+J  ARE  PLOTTED 
1,8H  AGAINST, 13, 31H  VALUES  OF  Ut.I  SCALE  FACTOR  *,E10.3] 

NP«NP 

DO  5  1=1,95 
5  PLOT { I J  sBLNK 
X  M  c  U  . 

X  M  U  *  0  . 

DO  10  1*1, NP 
1U  AlIUYlil) 

MSW*1 
GO  To  6  - 
20  DO  30  1=1, NP 
30  A  t I J  » Y3  t  I  ] 

M3W  =  2 
QO  TO  6  J 
40  DO  50  1=1, NP 
50  All  UY2S1  ) 

MSW  =  3 

60  DO  9 u  1=1, NP 

I F ( A  f 1 3  —  X M J  70,75.75 
70  IFfA[ I J-X MI )80,90,90 
90  XM  I  =  a  1 1  J 
GO  TO  90 
75  XM  =  A [  I  ] 

90  CONTINUE 

IFIMSW-IN-1) 320,40,100 
10-.  NRAS94-[N”1]*10 
RAsNRA 

SF* (XM-XMI ] /RA 

R6F»ABSF(XMI/[XM-XMin*RA*1.0 

N  R  r  R  F  F 
PRINT  3C- 

PRINT  31>- ,  NP ,  NP ,  SF 
GO  To(32>,  ,330,3401  ,N 
320  PRINT  325 

325  FORMAT  [  7X,  4MX  t .  J  ,  5X,  5HU  {  *  )] 

GO  TO  4 
330  PRINT  335 

335  FORMAT! 7X.4HX [ . ],5X,9HYlf*),5X,5HY?t=]  ] 

GO  TO  4 
340  PRINT  345 

345  FORMAT [7X,4HUt. ],2X,8HDY/DU[*]»4X,4HV[sJ,4X,4HW[*JJ 
4  PRINT  35 J 

35 U  F0RMAT(1X,116H . . . . 

1 . . . 

DO  290  I = 1 , N P 
PLOT  f NR] =PRD 


l»yiiI)/sf+ref 

PLOT ;L3  sSTAR 
GO  TO  1140,110,1101,9 
110  K*Y2f Il/SF+REF 
PLOT  f  K ) =£0 
IF [ N  — 2  3  150,150,120 
120  J»Y3f  n/SF  +  REF 
PLOTf JJsPLUS 
GO  Tn  l6d 

140  PRINT  7,Xtn,Yim,tPL0TtJl,J*l,95J 
PLOT ( L I *BLNK 
GO  TO  290 

150  PRINT  8.Xt!],Yltn,Y2Cn,  (PL0T(J1,J«1,85J 
PLOT  f LI *BLNK 
PLOT [ K ] =8LNK 
GO  TO  290 

16 u  PRINT  9,X[I],Yl[II,Y2tIJ»Y3t|l,[PL0TlJ],J*l,75] 
DO  165  N0=l,95 
165  PLOT ( NO ] =OLNK 
290  CONTINUE 
RETURN 
END 

02070,  17745 

PLOTT  0000012 
N  1000012 
NP  100001^ 

X  1000014  C  t  20144 
Y1  1000015  0020144 
Y2  1 C 0 0 0 16  0020144 
Y3  1000017  1)020144 
A  0 0 0  0036  0020144 
PLOT  0000346  0020144 
G  0017764 
BLNK  0017762 
PRD  0017760 
STAR  0017756 
EU  0017754 
PLUS  0017752 
COM  QU1775U 
EXE  0017746 
/ 0 0 0 U 7  05)00661 
/ 0  0  0  0 8  000067m 
/ 0  0  0  09  0000677 
/ 0 0 3 0 0  DC  3  0 70 6 
/  3  031 0  1)0  0  072 C 
I  000  r' 656 

+00001  00^0657 
-00096  0000775 
/000'j5  0001025 
XM  HfOlOOO 
0,  0001002 
XMU  0  0  0 1 0 -  4 
/ 0  0  0 1  0  001046 

-0  0  0  0  2  Q'J  01006 

MSW  0 0 0 1  10 
/00060  0OOH53 
/00 020  00" 10 74 
/00039  00 •.  10  76 
+  00002  0001  011 


/QQ34G  C0H124 
/9PG59  *001126 

*900 13  CO01312  1?2 

/0GC7C  C0-U167 
/00G75  0C?l23. 

X«!  3591014 
/9008C  0001231 
/90093  £(501216 
/30105s  C30l24t 

NnA  OCClOlA 
♦50C94  fio0lG17 
♦0001D  603102. 

RA  0001022 
SF  r-l  0125^ 

RcF  : 30126' 

48SF  l4i26 2  EXT  P3CG 

1.5  f i 0126* 

*R  000127- 
/■>03?0  06ei3*>5 
/’r33-:  D&61372 
/?034'S  OOCI^I'3 
/  ' 3325  0661361 
/3005,4  00^1437 
/  3335  C0fel376 
/'1345  C 0 ?  1 4 1 7 
/  -035C  56*1443 

L  -CO 1771 
/  ;  ?  1 4  p  CO  11632 
/SOUS  00H56. 

K  ..01273 
/'53155  ZC'Ji7j 
/3912€  00  161. 

J  .*01774 
/  *>916€  0601763 
/"C2q;  gG  ?C5f> 

-ft,*Cg6  C  C/ 1 275 
- 0 C C 7 6  ce  1276 
*.0  012  77 

/'C165  =  ?C<42 


SUBROUTINE  FNCT  (X,Y,F] 

COMMON  G>  dlNK> PRD » STAR* EQ> PLUS /COM*  EXE 

AArSORTF[2..1]  173 

8B*LOGF[G*EXPF(YJ-*1.0-G3  +  G*EXPFC-Y)*G 

CCaSoRTFtdB] 

F  *-A A*CC 

RETURN 

END 

00113,  17745 

FNCT  0000012 
X  1000012 
Y  1000013 
F  1000014 
G  0017764 
RLNK  OQ17762 

PRO  0017760 
STAR  0017756 

EQ  (i  6 1 7  7  5  4 
PLUS  0017752 

COM  ij  317750 
EXE  0  n  1 7  7  4  6 
A A  f 0  0  CO  '4 

SQRTF  0l4u02d  EXT  PROG 

2.C  000003* 

R3  C 000032 

LOGF  0143034  EXT  PROG 

EXPF  01 4*3  0 3o  EXT  PROG 

l.P  f 00004? 

CC  0000046 


oonooooci 


174 


C  **********  RUNG6-KU7TA  METHOD  FOR  SOLUTION  OF  O.D.E,  USING  GE225 
**********  H.  ALLEN  LINDSEY*  N0VEH8ER  1968 
**********  THE  FOLLOWING  VARIABLES  ARE  USED 
***  A  AND  B  ARE  THE  INTERVAL  LIMITS 

***  h  is  the  step  size 

***  y  IS  THE  INITIAL  VALUE  OF  THE  SOLUTION 
***  DELTA  IS  THE  INCREMENT  at  WHICH  the  SOLUTION 
IS  To  be^printed  out 
.**  g  is  the  Degree  of  ionization 

DIMENSION  NEQ 1 80  I # P  1 1 1 

DIMENSION  X1C200 J iYlf20Q I ,Y2t 2001 ,Y3 12001 
DIMENSION  Xll(70),Yllt7ol*Y22C705*Y33I70J 
COMMON. G,BLNK* PRD* ST AR, ?0, PLUS, COH, EXE 
999  FORMAT  isoaij 

998  FORMAT  I//5X*19h  ON  THE  INTERVAL  U**F8,4*6h  TO  U«*F8.4, 

1  17H  WITH  STEP  SIZE  »,F9.7/1 

997  FORMAT  [5X,51H  THE  TOTAL  NUMrER  OF  INTERVALS  FOR  THIS  PROBLEM  IS 
1  IA/1 

996  FORMiT  {5X,48H  THE  DEGREE  OF  IONIZaTION  FOR  THIS  PROBLEM  IS  S=* 

1  Fa. 4/1 

995  FORMAT  I5X,48h  ThE  EXPECTED  ERROR  POR  THIS  PROBLEM  IS  E*t+OR«)» 

1  El  6 . 6/ J 

994  FORMAT  18X,2H  U,l3X,2H  YS15X,2H  V//1 
993  FORMAT  I5X,F12.7*3X,F13.31 
992  FORMAT  [ 5X, Fi2 . 7 , 3X , Fl3 , 8 > 3X, Fl3 . 8 1 
<391  FORMAT  [6F10.3] 

99C  FORMAT  I//////1 

989  FORMAT  I2X,1A1,2X*1A1*2X*1A1*2X*1A1*2X91A1] 

988  FORMAT  [ 57 , 37H  U*NORMALjZED  DISTANCE  FROM  JUNCTION  / 

1  5X* 33H  YzNORmALIZED  ELECTRON  POTENTIAL  / 

2  5X,35H  VzNORMALlZED  SPACE  CHARGE  DENSITY  / 

3  5X*32H  WzSCHOTTKY  APPROXIMATION  FOR  V  ////} 

READ  989>PRD*STaR*Eq,PLUS*BLNK 
Nl  =  3 

C  **********  MAI<^  PROGRAM 

1  READ  991*A,8,DFLTA,H,Y*G 
f  **********  CHECK  FOR  EOT 

IF  { H J  2,100,2 

2  NHrtR-AJ/H 
PRINT  99o 

READ  999,  INEOin  ,  1=1,80] 

E*w**5 . 0 

P[1}=SGRTF{2,0*Y1 
c  **********  print  headings 
PRINT  999  *  NEO 
PRINT  998»A,B,H 
PRINT  997 * NH 
PRINT  996* G 
PRINT  995 *  E 
PRINT  988 
PRINT  994 
PRINT  993, A, Y 

C  **********  INITIALIZATION 
10  X  *  A 
1  =  1 

J* 

C  **********  COMPUTATION  OF  THE  NUMBER  OF  PRINT  INTERVALS 
N*DEl.  TA/H 


j-o 


nn*n*i 
15  Y0*Y 

»**•****«»  EVALUATION  of  The  function  for  runge-kutta  forhula 

30  CALL  FNCT  (X,YO,F) 

ZK1  =  F 

XH=X*H/2 . C 
Y0=Y*H*ZKl/2. 0 
CALL  FNCT  [XH,Y0/FJ 
ZK2-F 

YO=Y*H*ZK2/2.0 

CALL  FNCT  IXH,YO,FJ 

ZK3*F 

XHr  X-t-H 

Y0=Y*H*ZK3 

CALL  FNCT  [XH,YO,F] 

ZK4s  f 

50  YsY*M*[ZKl+2,0*lZK2*ZK3l*ZK4J/6«0 
D  =  l.''*I1.0/G-1.01*EXPFt-Y] 

V*l. J/0-0*EXPFC-Y) 

**********  INCREMENT  ARGUMENT  AND  COUNTERS 
7*1 

X  *  A  7  *  H 
1*1*1 

***(******  TEST  FOR  PRINT  INCREMENT 
IF  I  I -UN J  30,70,70 
7 1‘  PRINT  992,X,Y,V 
J  =  J-1 

Y1?J1=Y 

Y2 ( J 1  - V» Y1 { 11 

IF  IX-PllJ}  76,76/78 
76  Y3iJl=YlU) 

CO  TO  79 

78  Y3fJj=C,. 

79  X1(J]=X 

**********  TEST  FOR  FINAL  COMPUTATION 
IF  [mN- (NH+lS 1  S0.9C.90 
8"  NNsN+NN 
GO  TO  15 

**********  STORAGE  /  SCALE  FOR  PLOT  SUBROUTINE 
90  Ss  1 1 . 0/DELT  A  3 
MN  =  S 

IF  IS-58.3J  96,96/92 
9?  SCALF=S/58. f*i. o 

nm=scale 

L  =  m.\/nM 
K=-2*NM 

DO  94  JJ  =  1 , L  / 1 
KsK+NM 

YllljJ)=Yl{K+NM*lJ 
Y2?lJJ33Y2iK^NM+iJ 
Y33 [ Jj ] =Y3 { K+NM+1 J 
94  Xll{jJI=Xli«+NM+l] 

NP  =  L 

GO  TO  9‘3 

96  CONTINUE 

DO  97  J1=1,MN,1 
Yll(Jl)=Yli J1I 
Y2?i  J1]=Y2{ J1 ] 

Y33 [ Jl I =Y3  C  J1 } 

97  Xll  l  J1J=X1U1J 


nooo  ooooo 
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**********  RUNGE-KUTTA  METHOD  'OR  SOLUTION  OF  O.D.E.  USING  GE225 
**********  H,  ALLEN  LINDSEY>  N0vEM8cR  1968 
**********  THE  FOLLOHING  VARIABLES  ARE  USED 
***  A  AND  B  ARE  THE  INTERVAL  LIMITS 

***  h  is  the  step  size 

***  Y  IS  THE  INITIAL  VALUE  OF  THE  SOLUTION 
***  DELTA  IS  THg  INCREMENT  AT  WHICH  THE  SOLUTION 
IS  TO  reprinted  OUT 
***  g  is  the  degree  of  ionization 

DIMENSION  NEQ  [ 80  $  «  P ( 1 1 

DIMENSION  XlI200  UYlt23cUY2t200  3  >Y3I2003 
DIMENSION  Xli{70J*Yiit70l,Y22C70J»Y33(70J 
COMMON  G,8LNK,PRD»STAR«g0, PLUS, COM, EXE 
999  FORMAT  1 8 0 A 1 1 

998  FORMAT  1//5X,19h  ON  THE  INTERVAL  U*»F8.4,6h  TO  Ui,F8.4, 

1  17H  HUH  STEP  SIZE  a,F9.7/5 

997  FORMAT  l5X,51H  THE  TOTAL  NUmrER  OF  INTERVALS  FOR  THIS  PROBLEM  IS 
1  16/ 1 

996  FORMaT  {5X,48H  THE  DEGREE  OF  IONIZATION  FOR  THIS  PROBLEM  IS  G=, 

1  Ffl.4/1 

995  FORMAT  {5X,48h  TkE  EXPECTED  ERROR  FOR  THIS  PROBLEM  1$  E«UOR*U 
1  E j  6 . 8/ ] 

994  FORMAT  I8X,2H  U,i3X,2H  Y,l5X,2H  V//) 

993  FORMAT  I5X,Fl2.7,3X,rl3,9l 

992  FOPMaT  [5X,Fi2.7,3X,ri3,8,3X,Fl3.83 

091  FORMAT  I6F10.3] 

990  FORMAT  I//////] 

989  FORMAT  {2X,1A1,2X,1A1»2X,1A1,2X,1A1,2X»1A1] 

988  FORMAT  {5X,37H  U*N0RMALIZED  DISTANCE  FROM  JUNCTION  / 

1  5X,33H  Y*N0RMALIZED  ELECTRON  POTENTIAL  / 

2  5X » 35H  V*N0R.MALI2ED  SPACE  CHARSE  DENSITY  / 

3  5X,32H  WzSCHOTTKY  APPROXIMATION  FOR  V  ////J 
READ  989, PRD, STAR, Eq,PLUS,BLNK 
Nl  =  3 

C  **********  MAIN  PROGRAM 

1  READ  991,A,8,DFLTA,H,Y*G 
C  **********  CHECK  FOR  EOF 

IF  £  M )  2,100,2 

2  NH*IB-A]/H 
PRINT  99o 

READ  999, {NEGCIJ, 1=1,80] 

P*W**5. 0 

P(i)=SQRTF[2.0*YJ 

C  **********  print  headings 

PRINT  999, NEQ 
PRINT  998,A,B,H 
PRINT  997, NH 
PRINT  996, G 
PRINT  995, E 
PRINT  988 
PRINT  99 4 
PRINT  993, A, Y 

C  **********  INITIALIZATION 
10  X  =  A 
1  =  1 
J* 

C  **********  COMPUTATION  OF  ?WE  NUMBER  OF  PRINT  INTERVALS 

n*dei.ta/h 
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NN=N*1 
15  YO*Y 

c  **********  EVALUATION  or  the  function  for  RUNGE-KUTTA  formula 

30  CALL  FNCT  [X,YO,FJ 
ZK1  =  F 

XH=X*H/2. C 
YO=Y+H*ZK1/2.0 
CALL  FNCT  (XH,YO»FJ 
ZK2  =  r 

Y0sY*H*ZK2/2, 0 
CALL  FNCT  tXH,YQ,F] 

ZK3*F 

XH=X*H 

Y0=Y*K*ZK3 

CALL  FNCT  [XH,YO,FJ 

ZK4  =  r 

5 0  Y=Y*M*[ZKl+2,0*tZK2*ZK3]iZK4]/6»0 
D  =  1.0-'(1.0/G-1.03*EXPFt-yJ 
V*l.  )/D-G*EXPFt-YJ 

C  **********  I NCREMENT  ARGUMENT  AND  COUNTERS 
Z*I 

X* A+7*H 
I  *  I  ♦! 

c  ***;*•**•«  TEST  for  PRINT  INCREHENT 
IF  c  T -M-M 3  30,70,70 
7,.  PRINT  992 ,  X,  Y,  V 
J  =  J*1 
Ylf J}=Y 
Y2  f  J 1 =  V*Y1 ( 1 J 
IF  IV-PIlJ]  76,76/78 
76  Y3 I J ] - Y1 [ 1 ] 

GO  TO  79 

78  Y3fJj=G,*; 

79  X1IJ]=X 

C  **********  TEST  FOR  FINAL  COMPUTATION 
rF  IMN-[NH*1J]  30,90,90 
8 NNsN*NN 
GO  To  15 

C  **********  STORAGE  /  SCALE  FOR  PLOT  SUBROUTINE 
90  Ss 1 1 . 0/DELT  A  3 
MN  =  S 

IF  IS-58.DJ  96,96,92 
9?  SCALP=S/58.n*i.o 
NM=SCALE 
L=MN/NM 
k=-2*NM 

DO  94  JJ=1,L> 1 
K*K+NM 

Ylll  JJJsYltK*NM  +  lJ 
Y22I JJ]3Y2(k*NH*1] 

Y33(JJ]=Y3{K<-NM*1J 
94  XU  [  JJ)=Xl[K  +  ivM  +  i] 

NP  =  L 

GO  TO  98 

96  CONI  [Nile 

OO  97  J1=1,MN,1 
Yllf JlJsYllJlJ 
Y22I Jl)=Y2t J13 
Y33 I Jl ] = Y3 ( J1  ] 

97  XIII Jl] =X1{ Jl] 
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np=m.*j 

9«  CALL  DLOT  T  (  Ml ,  NP>  XU ,  Yli,  Y22 >  Y33  ] 
GO  Tn  1 
100  CALL  EXIT 

pno 

■‘5616,  17745 

,KP  coo-  •  :■ ) 

MEQ  000  0.1  0020120 

P  000  124  0J20001 

XI  000  126  0020310 
Y1  000  746  0520110 
Y2  0001566  0 2 0 3 1 G 
Y3  00024:^6  '.’^20310 

XII  0003226  00201., 6 
Yll  O0M344?  0.:291T6 
Y22  0  3  03636  U.S20U6 
Y33  n004>7?  0 r< 2 0 1 C 6 
G  0017764 


PLNK 

17762 

PRO  0 C 1 7 7 6 0 

STAR 

17756 

EO  U  il 

7754 

PLUS 

17752 

CO'H  i  017751 

EXE  r.<ii7746 

/  -09Q9 

r  .4321 

/  •’  0  9  o  8 

4324 

/00997 

•  i.  4356 

/00996 

•  5.4404 

/  1Q9Q5 

•  f  4  4  3 1 

/ '!  0  90  4 

’  445  7 

/•'  G9Q3 

C>4471 

/  -3902 

0451  ) 

/  o  n  o  9 1 

•4512 

/  10O9Q 

T4515 

/ 00989 

u45?  : 

/)|09R8 

453? 

Mi  f.  0  0  4  3  r.  ^ 

•‘-00003 

■  <14307 

/10001 

j4644 

A  0  0  0411)' 

R  000431? 

DELTA  ,4314 

H  0  r  0  4  3 1 4 

Y  -..(0466 
/ 1 0  0  0  2  ->472 

/"01C0  -1.5614 

■’:H  001466? 

i  onn46^A 

- r  r  - 1  4^47 

-■’CCPl  -!467 
P  001467? 

5.0  i, 0  04674 

,19  (.1^4474  PIT  py.jR 
SQRTP  14470.-  EXT  D90G 

2.0  ''0  047'? 

.AY  i-1047.)4  EXT  PROG 
/ r  C  0 1 0  05032 

X  .115  05  3  4 


J  U'J05  >36 
♦  'JOOOO  05037 
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FNCT  1  i'S044  EXT  PROG 

F  C  0  0 5  >46 
ZK1  0005^50 
XH  0005-5? 
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ZK4  i)005  6r 
/ 0  0  0  5  0  05165 

6.0  0005201 
D  u "05202 

1.0  00052j4 

E XPF  145206  EXT  PROG 

V  000521? 

Z  0005214 


/00070 

i5304 

/00076 

0533/ 

/O0078 

(15346 

/90079 

15354 

0.0  000522n 

/30080 

i  5373 

/D0090 

j5377 

S  000522 2 

MN  u  0  05224 

58.0 

u  5  4 1  [i 

/30G96 

?  5  5  5 1 

/0  0  092 

■  o544o 

SCALE 

.,-'■5412 

NM  0005414 

L  0005415 
K  0005416 
-  0  0  0  U  2  ■>541/ 

JJ  .j 0  054 2  i 
/00094  1 V  5  5  2  4 

NP  0005422 

/  1  n  0 g 8  '  56 G 4 

J1  0005421 
/0Q097  15565 

PLOTT  1"54?4  EXT  PROG 

EXIT  1  ,542o  FXT  PROG 


nN  me  interval  u*  n.oooo  to  u*  20,0000  with  step  size  «o. 200000 

THE  total  NUMdFR  OF  INTERVALS  FOR  THIS  PROBLEM  IS  iOO 

the  degree  qf  Ionization  for  this  problem  is  g*  0.0500 

THE  EXPECTED  ERROR  for  THIS  PROBLEM  IS  E*!>OR-]  0.31999996E-03 

UsNORMALlZED  DISTANCE  FROM  JUNCTION 
YrNORMALlZED  ELECTRON  POTENTIAL 
V*NfiRMAL  I  ZED  SPACE  CHARGE  DENSITY 
WsSCwOTTKY  approximation  FOR  V 
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P.POO  1  1G0 
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9,?000./ 00 
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9.4000000 

3,75659978 
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1  '.2000000 

2.77593549 

1  ,4OOUu00 
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2.24085804 

11 ,  y00  .’"00 

2.09067604 

11  ,?00'Jj00 

1.95220643 

11,4000  0  O'  0 

1.82428385 

11,6000 w 00 

1.70588840 

ll.fiOUJoOO 

1.59612629 

12.000 ) ! 00 

1.49421241 

12,20  0  j'.  0  0 

1.39945500 

12.400  :•  .00 
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13,600"' t/OO 

0.8927Q322 
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0.837898/3 

14  ,  OOO'J.jOG 

0.78656939 
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0.69339474 
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15.n00ji.0  0 
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16,'jGO  .00 

0.41995577 
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16,600v  0  0 

0.34818255 

16,80  0  t o 0 0 

0.32710756 

16,9999999 

0.30731226 

17,?00v  < 0 0 

0.20871801 
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0.27125119 

17,600')/  0-J 

0 .25484284 

1 7 , 8  0  n  i  j  0  1 

0.23942836 

18,00  0  JoOO 

0.22494717 

18,20  0  ’00 C 

0.21134246 

18,3999999 

0 .19856097 

18, 600.1  J 00 

0.18655271 

1 8 , 0  0  0  I,  j  0  0 

0.17527077 

13,999999° 

0.16467112 

19,200  1000 

0.15471242 

19.40  0  /  -  0  0 

0.14535586 

1° . 6  0  0  /  1 0  0 

0.13656501 

19.80  0  .  .00 

0.12830563 

19,9999999 

0.12054558 

0.81068663 

0.75317845 

0.69140161 

0.62833644 
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0.50842300  180 

0 . 45485896 

0.40655608 

0.36358749 

0.32570041 

0.2924639? 

0.26337331 

0.23791585 

0.21560783 

0.19601236 

0.17874544 

0.16347535 

0.14991866 

0.13783483 

0.12702U43 

0.11730378 

0.10853994 

0.10060648 

0.09339975 

0.08683184 

0.08082791 

0.07532408 

0.07026562 

0.06560544 

0.06130290 

0.05732273 

0.05363422 

0.05021051 

0.04702801 

0 .04406591 

0.04130575 

0.03873111 

0.03632733 

0.03408124 

0.03198099 

0.03001585 

0.02817610 

0.02645290 

0.02483815 

0,02332445 

0.02190502 

0.02057358 

0.01932438 

0.01815206 

0.01705170 

0.01601869 

0.01504879 

0.01413802 

0.01328268 

0.01247934 

0.01172476 
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A  theoretical  study  of  metal-semiconductor  (It-type)  rectifying  contacts 
is  developed.  This  study  begins  by  first  analyzing  previous  models  for 
this  type  of  junction.  Particular  attention  is  given  to  the  Scbottky 
model  and  to  the  approximations  it  contains.  This  node!  is  then  improved 
upon  by  taking  into  account  nonuni form  impurity  ionization  and  the  free 
electron  concentration  in  the  depletion  region.  Using  this  more  exact 
model  a  theoretical  expression  for  the  differential  junction  capacitance  is 
calculated.  The  results  indicate  that  the  junction  capacitance  as  a  function 
of  reverse  bias  can  be  used  to  accurately  predict  the  doping  concentration  in 
the  semiconductor  material,  but  does  not  yield  a  correct  measurement  of  the 
equilibrium  diffusion  potential  or  barrier  height. 

The  current  voltage  characteristic  for  this  type  of  contact  is  also  discussed 
An  expression  for  the  I-V  characteristic  of  this  junction  is  derived  based 
upon  a  diffusion  model.  This  expression  is  then  improved  upon  by  accounting 
for  tunneling  and  quanfcum-mcchanical  reflection  of  carriers  at  the  junction. 


